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Abstract 



In one-dimensional quantum wires the interplay of electron correlations and impu- 
rities strongly influences the low-energy physics. The diversity of energy scales and 
the competition of correlations in interacting Fermi systems can be treated very 
efficiently with the functional renormalization group (fRG), describing the gradual 
evolution from a microscopic model Hamiltonian to the effective low-energy action 
as a function of a continuously decreasing energy cutoff. The fRG provides the uni- 
versal low-energy asymptotics as well as nonuniversal properties, and in particular 
an answer to the important question at what scale the ultimate asymptotics sets in. 

The lowest order truncation of the fRG hierarchy of flow equations consid- 
ered previously for spinless fermions is generalized to spin-| systems and extended 
including renormalization of the two-particle interaction, in addition to renorma- 
lization of the impurity potential. The underlying approximations are devised for 
weak interactions and arbitrary impurity strengths. A comparison with numerical 
density-matrix renormalization results for systems with up to 1000 sites shows that 
the fRG is remarkably accurate even for intermediate interaction strengths. 

We investigate the influence of impurities on spectral and transport properties 
of fermionic lattice models with short-range interactions. The results capture rele- 
vant energy scales and crossover phenomena, in addition to the universal low-energy 
asymptotics. For weak and intermediate impurity strengths the asymptotic behav- 
ior is approached only at rather low energy scales, accessible only for very large 
systems. For spin-i systems two-particle backscattering leads to striking effects, 
which are not captured if the bulk system is approximated by its low-energy flxed 
point, the Luttinger model. In particular, the expected decrease of spectral weight 
near the impurity and of the conductance at low energy scales is often preceded by 
a pronounced increase, and the asymptotic power laws are modifled by logarithmic 
corrections. 
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1 Introduction 



In one dimension metallic electron systems are strongly affected by interactions. Dif- 
ferently from the conventional Fermi-liquid behavior, t he generic low-ene rgy physics 



is described by the Luttinger- liquid phenomenology [Giamarchil |2004| . For vari- 
ous correlation functions Luttinger-liquid theory predicts anomalous power laws; 
for spin-rotation invariant systems the exponents can be expressed in terms of a 
single interaction-dependent parameter Kp. An important aspect concerns the pe- 
culiar effects due to the interplay of impurities and interactions. For Luttinger liq- 



uids with repulsive interactions a . 



Luther a,nd Peschellll 974 iMatti: 



ready a single static impurity has a strong effect 



19741 : lApel and R.icelll982l : iGiamarchi and Schulz 



Il988| . At low energy scales even a weak impurity effectively "cuts" the system into 
two parts with open boundary conditions at the en d points, and physical o bserv- 



ables are controlled by the open chain fixed point [Kane and Fisher! Il992al jcl|. In 
particular, the impurity potential becomes dressed by long-range oscillations lead- 
ing to a characteristic power-law suppression of the local density of states near the 
impurity and the conductance through the impurity down to zero in the low-energy 
limit. The asymptotic behavior is universal in the sense that the exponents depend 
only on the properties of the bulk system via Kp, while they do not depend on the 
impurity strength or shape. These power laws are generally modified by logarithmic 
corrections in the presence of two-particle backscattering. 

The asymptotic low-energy properties of Luttinger liquids with a single impurity 
are rather well understood. Universal power laws and scaling functions have been 
obtained by bosonization, conformal field theo rv and exact sol utions for the low- 
Numerical meth- 



energy asymptotics in special integrable cases 
ods as exact diagonalization and the density-matrix renormalization group (DMRG) 
confirm the field-theoretical predictions and the validity of the underlying assump- 
tions for microscopic fermionic systems with Luttinger-liquid behavior. The limited 
system size accessible to numerical solutions is however a serious constraint for a 
systematic analysis beyond the perturbatively accessible weak and strong-impurity 
regimes. The important question arises at what scale the ultimate asymptotics sets 
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1 Introduction 



in and asymptotic power laws are actually valid. That scale can indeed be surpris- 
ingly low, and the properties above it very different from the asymptotic behavior. 
Recently a functional renormalization group (fRG) method has been introduced for 
a direct treatment of microscopic models of interacting fermions, which does not 
only capture correctly the universal low-energy asymptotics, but allows to compute 
observables on all energy scales, providing thus also nonuniversal properties, and a 
possible key to the understanding of the behavior at intermediate scales accessible in 
experiments. Some of the nonuniversal properties can be computed numerically by 
the DMRG, but this method is limited to lattice systems with about 1000 sites, and 
only a restricted set of observables can be evaluated with affordable computational 
effort. 

The fRG provides a powerful computational tool to study interacting Fermi 
systems, especially low-dimensional systems with competing instabilities and en- 
tangled infrared singularities. Starting point is an exact hierarchy of differential 
flow equations for the Green or vertex functions of the system, describing the grad- 
ual evolution from the microscopic model Hamiltonian to the effective action as 
a funct ion of a continuo usly decreasing energy cutoff introduced in the free prop- 
agator Salmhoferl Il998l | . Approximations are then constructed by truncating the 



hierarchy and parametrizing the vertex functions with a manageable set of variables 
or functions. The fRG captures the expected universal power laws at low energy, 
as well as relevant energy scales and nonuniversal crossover phenomena at inter- 
mediate scales, as for the temperature dependence of the conductance through a 
double barrier 



Enss et al 



2005; 



Meden et al. 



200q . The direct application to mi- 



croscopic models allows for a flexib 
rings threaded by a magnetic flu x 



Barnabe-Theriault et al 



2005 



3b| 



e modeling of different geomet ries, as mesoscopic 
Meden and SchoUwock 2003a |b| and Y junctions 



Inprevious applications to spinless Luttinger liquids with impurities 



Meden et al, 



the fRG hierarchy of flow equations was truncated at first order, where the 
renormalized vertex is approximated by the bare interaction. Despite the simplicity 
of this scheme the effects of a single static impurity are captured qualitatively, and 
for spinless fermions in the weak coupling limit also quantitatively. It turned out 
that the asymptotic behavior typically holds only at very low energy scales and for 
very large systems, except for very strong bare impurities. 

In the present work we further develop and extend the fRG approach for Lut- 
tinger liquids with impurities to spin-i fermions and include two-particle vertex 
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renormalization, in addition to the renormalization of the impurity potential The 
underlying approximations are devised for weak interactions and arbitrary impu- 
rity strength. A comparison with exact numerical DMRG results for systems with 
up to 1000 sites shows however that the fRG with the inclusion of vertex renor- 
malization is remarkably accurate even for intermediate interaction strengths. For 
spinless fermions this extension improves considerably the quantitative accuracy of 
the results in particular at intermediate interaction strengths, whereas for spin-| 
systems vertex renormalization is necessary to take into account that backscatter- 
ing of particles with opposite spins at opposite Fermi points scales to zero in the 
low-energy limit. Explicit flow equations are derived for various lattice fermion 
models supplemented by different types of impurity potentials. We present results 
for spectral properties of single-particle excitations, the oscillations in the density 
profile induced by impurities or boundaries and the linear conductance for chains 
with up to 10^ lattice sites. Two-particle backscattering leads to peculiar effects, 
which are not captured if the bulk system is approximated by its low-energy fixed 
point, the Luttinger model. In particular, the expected decrease of spectral weight 
near the impurity and of the conductance at low energy scales is often preceded by 
a pronounced increase, and the asymptotic power laws are modified by logarithmic 
corrections. 

The outline of the thesis is as follows. 



In Chapter |21 we give a short overview on general aspects of Luttinger liquids 
with impurities. 

The fRG formalism is developed in Chapter El We briefly review the fRG 
for interacting Fermi systems, and derive the hierarchy of differential flow 
equations for the one-particle irreducible (IPI) vertex functions. 

In Chapter E] we describe the implementation of the fRG technique for vari- 
ous one-dimensional microscopic lattice models with impurities, providing de- 
tails on the parametrization of the two-particle vertex and different truncation 
schemes. Parts of this chapter are published in 

S. Andergassen, T. Enss, V. Meden, W. Metzner, U. Schollwock, and K. Schon- 
hammer. Functional renormalization group for Luttinger liquids with impuri- 
ties, |Phys. Rev. B 70, 075102 (2004)| |cond-mat/04035171 
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1 Introduction 



S. Andergassen, T. Enss, V. Meden, W. Metzner, U. Schollwock, and K. Schon- 
hammer, Renormalization group analysis of the one- dimensional extended Hub- 
hard model with a single impurity, |cond-mat/ 05090211 

In Chapter we present results for spectral properties of single-particle ex- 
citations near an impurity or boundary, the density profile, and transport 
properties in the presence of a single and a double barrier. In the first part we 
focus on spinless fermions; the modifications due to the spin degree of freedom 
are addressed in the second part. Parts of this chapter are presented in the 
above publications, and for the conductance in 

V. Meden, S. Andergassen, W. Metzner, U. Schollwock, and K. Schonhammer, 



Scaling of the conductance in a quantum wire, Europhys. Lett. 64, 769 (2003 
lcond-mat/ 0303460'; 

V. Meden, T. Enss, S. Andergassen, W. Metzner, and K. Schonhammer, 
Correlation effects on resonant tunneling in one- dimensional quantum wires, 
jPhys. Rev. B 71, 041302(R) (2005)1 lcond-mat/ 0403655l 

T. Enss, V. Meden, S. Andergassen, X. Barnabe-Theriault, W. Metzner, and 
K. Schonhammer, Impurity and correlation effects on transport in one-dimen- 
sional quantum wires, [Phys. Rev. B 71, 155401 ( 2005)| |cond-mat/0411310t 



S. Andergassen, T. Enss, and V. Meden, Kondo physics in transport through 
a quantum dot with Luttinger liquid leads, .,cond-mat/0509576i 

We conclude in Chapter IHlwith a summary and an outlook on further appli- 
cations and extensions of the present work. 
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2 Impurities in Luttinger liquids 



The exactly soluble Luttinger model provides a generic scenario for one- dimensional 
Fermi systems with repulsive interactions, denoted as "Luttinger liquid". The low- 
energy physics is completely determined by a few interaction- dependent characteristic 
parameters describing the power-law exponents of the correlation functions. Already 
a single static impurity leads to peculiar modifications of the electronic properties of 
Luttinger liquids. Even for a weak impurity potential, physical observables behave 
as if the system is split into two parts in the low- energy limit. The local density of 
states near the impurity and the conductance through the impurity vanish as power 
laws. 



2.1 Luttinger liquids 



In one-dimensional interacting Fermi systems Fermi-liquid theory is not valid. The 
breakdown of Fermi-liquid theory is indicated already in second order perturba- 
tion theory, where the reduction of the quasi-particle weight at the Fermi surface 
due to interactions diverges logarithmically. These divergencies can be treated by 
a weak-coupling renormalization- group raethod applied to an effective low-energy 
theory known as g-o\ogy model [Solyom I1974 Depending on the values of the 
bare couplings, the renormalized couplings flow either to strong coupling, and hence 
out of the perturbatively control 
exactly soluble Luttinger model 



led regime, or to a fixed-point Hamiltonian, the 


Tomonaga 


mi 


Luttinfferl 




Mattis and Lieb 



The term "Luttinger liquid" has been introduced for the latter systems, in 



analogy with the mapping of low-energy states of inte racting e 
the Fermi gas in higher dimensions for Fermi liquids Haldane 



ectron systems onto 



1981b|. The normal 
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2 Impurities in Luttinger liquids 



gapless metallic phase is characterized hj i) a. continuous momentum distribution 
with a power-law singularity at the Fermi surface, described by a nonuniversal ex- 
ponent a; ii) a single-particle density of states which vanishes as near the Fermi 
energy, implying the absence of fermionic quasi-particles; Hi) finite charge and spin- 
density responses for long wavelengths and the existence of collective bosonic charge 
and spin-density modes; iv) power-law singularities in various correlation functions 
with interaction-dependent exponents; v) separation of spin and charge degrees of 
freedom. There are several go od reviews on one-dirnension al Fermi systems, recent 
reviews are presented in Refs. Voitlll995c lGiamarchill2004| . In the following we will 



summarize the most important results. 

Theoretical work on interacting fermions in one dimension has progressed along 
different lines . Besides the perturbative investigation of the weak-coupling limit 



S61voJll979t . Luttinger-liquid theory is usually formu lat ed using the bosonization 



technique Mattis and Lieb 



tecnnique |M 
lMattijh974( . 



integrable models 



1965; 



Haldane 



1980, 



198l3Lt 



Luther and Peschel 



1974 : 



A dif ferent approach i s based on the Bethe-ansatz method for special 



The computation of correlation functions is 



however very difficult from the complicated expre ssions for the eig enfunctions. 



Haldane 



Il981b| . 



the low-energy 



As proposed in a seminal work by Haldane 
physics of the Luttinger model is generic for interacting fermions in one dimen- 
sion with repulsive interactions. In the language of the renormalization group the 
Luttinger model Hamiltonian is the fixed-point Hamiltonian for a large class of one- 
dimensional fermions with repulsive interactions. The Luttinger model can be solved 
exactly at any interaction; it is characterized by a linear dispersio n relation, and th e 
electron-electron interaction is limited to forward scattering only |Giamarcliill2nn4| . 
Umklapp and backscattering processes, as well as additional terms for more general 
models arising f rom band curvature are irrelevant and vanish in the low-energy limit 



Haldane 



1981bl | . As in the Landau Fermi liquid a few parameters completely deter- 



mine the low-energy physics. The charge degrees of freedom of Luttinger liquids are 
described by a sound velocity Vp and the dimensionless parameter Kp, and the spin 
degrees of freedom are characterized by a spin- wave velocity v^j and K^. All corre- 
lation functions are uniquely parametrized by K^, and the velocities of the collective 
modes f ^, with u = p,a, in the low-energy limit; the corresponding exponents are 
determined by K,^. For noninteracting particles Kp = K„ = 1. In the absence of 
a magnetic field, the ground state is spin-rotationally invariant and = 1, while 
Kp < 1 (> 1) for repulsive (attractive) forces. 
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2.1 Luttinger liquids 



For = 1 the momentum distribution function exhibits a power-law singular- 
ity at the Fermi level with e xponent a = {K^ + — 2)/4 for any nonvanishing 



interaction IGiamarchill2C 

kp obeys a power law 



For a < 1 the momentum distribution function near 



\n{k)-n{kF)\ ~ \k - kpl"" . 
The spectral function has the form 



(2.1) 



(2.2) 



in the low-energy limit. Landau quasi-particle excitations are absent in the Luttinger 



liquid. The power la ws hold a 



nonlinear dispersion Haldane 



so for nonsoluble generalizations of the model with a 



1981b | 



In the Luttinger model the charge and spin density modes are exact undamped 
eigenstates, and any excited state of the model is a superposition of these elementary 
excitat ions. This becomes par ticularly explicit in the bosonized form of the Luttinger 
model Mattis and Lieblll965l |. The Luttinger model Hamiltonian conserves charge 



and the z component of spin separately on each Fermi point. Charge and spin 
excitations are completely independent, as the respective terms i n the Hamiltonian 
commute. This phenomenon is called "spin-charge separation" iGiamarchilbnn^ . 
charge and spin propagate with different velocities. 

Concerning the leading low-energy long-wavelen gth response func tions there is 
no difference between Fermi and Luttinger liquids 



Giamarchi 



200J. Thermody- 



namic properties as the compressibility and the susceptibility do not differ from 
the Fermi-liquid description and the modification due to the interaction leads to 
renormalized coefficients depending on and v^. Differences between Fermi and 
Luttinger-liquid behavior arise only from the enhanced phase space for forward scat- 
tering in one dimension. Marked differences appear in the single-particle propagator, 
which determines the momentum distribution function and the spectral density for 
single-particle excitations. In a Fermi liquid residual interactions modify the propa- 
gator only on a subleading level, leading for example to a small quasi-particle decay 
rate, while in a Luttinger liquid forward scattering affects the leading low-energy 
behavior. Another distinctive feature of Luttin ger liquids is the singular behavior 
of density correlations with momenta near 2kF loiamarchilEc ' 



Conservation laws play a crucial role in one-dimensional Fermi systems 



Metzner et al. 



|l998[. In addition to the usual charge and spin conservation the discrete structure 
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2 Impurities in Luttinger liquids 



of the Fermi surface in one dimension leads to an additional conservation law: sep- 
arate charge conservation in low-energy scattering processes for particles near the 
left and right Fermi points, respectively. Separate spin conservation is spoiled by 
the backscattering process generally present in models of spin-| fermions. In most 
cases of interest, in particular for the models considered in the present work, the 
backscattering amplitude scales to zero at low energies, and the separate spin conser- 
vation is restored asymptotically. The velocities associated with the corresponding 
conserved curre i its provide a complete param etrization of the low-energy physics 



Haldane 



1981b: 



Metzner and Di Castro 



1993 1 . 



2.2 Impurity effects 



An important aspect of Luttinger-liquid behavior concerns the peculiar modifica- 
tion of the electronic properties in the presence of impurities. For Luttinger liq- 
uids with repulsive interactions {Kp < 1) already a single static impurity has a 



Luther and Peschel 


1974; 


Mattis 


1974; 


Auel and Ric€ 


1982: 


Giamarchi and Schulz 


1984 Kane and Fisher! 


1 992alf : iFurusaki and Naeaosa 


1993a 


b: Yue et al .111 9941. In 



general the interplay of disorder and interactions is still a challenging issue, although 
the properties of noninteracting disordered electronic systems are rather well under- 
stood. In one-dimensional noninteracting systems disorder leads to localization of 
all electrons; the localization length characterizing the exponential deca y of the 
wave function is of the same order as the mean free path |Giamarcliill2nn4| . On the 
other hand interactions strongly affect the properties of the pure system, leading 
to Luttinger-liquid behavior. Thus in one dimension a particularly strong mutual 
influence of disorder and interactions is expected. 

Relevant parameters in the description of disorder are the strength of the indi- 
vidual impurity V and the impurity density nimp- The variation of these two param- 
eters leads to different physical effects. In the limiting case of very weak individual 
impurities with a dense distribution the effect of a single impurity is negligible and 
collective effects dominate; the corresponding relevant length scale is ~ 1/nimp- As 
a consequence of the central-limit theorem, for continuous systems the disorder can 



be described by a Gaussian distribution in the limit n 



stant rZimpV^^ measuring the disorder strength (Giamarchil |2004| . The main results 



oo and — ^ for con- 
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2.2 Impurity effects 



: or Gaussian diso rder from a perturbative treatment can be summarized as follows 
Giamarchil 120041 . Interactions are effectively renormalized by disorder, which is 
reversely affected by interactions. Repulsive interactions generally enhance local- 
ization whereas attractive ones reduce this effect. For spinless fermions attractive 
interactions enhance superconducting fluctuations, leading to an effective screen- 
ing of the disorder. For spin-| fermions a competing effect arises. The tendency 
towards a uniform charge distribution inhibits the coupling to disorder, leading to 
an increase in the localization length for strong interactions in the pure Hubbard 
model. For the extended Hubbard model with a local as well as nearest-neighbor 
interaction this effect is reduced. The opposite limit examined in the present work 
corresponds to strong and dilute impurities. In this case collective effects do not 
play any role and the problem essentially reduces to a single isolated impurity. An 
interesting unsolved problem concerns the combination of single impurity and col- 
lective effects at intermediate scales: depending on whether collective effects become 
important before the individual impurities renormalize to high barriers, a different 
characteristic behavior is expected. 

In the following we consider the case of a single or double impurity, where the 
effects of coherent scattering from many impurities are absent. The asymptotic 
low-energy properties of Luttinger liquids with a single impurity have been inves- 
tigated by mapping the problem onto an effective field theory, where terms which 
are expected to be irrelevant in the low-energy limit are neglected. For attractive 
interactions the impurity is irrelevant in the renormalization-group sense and scales 
to zero at low energies. For repulsive electron systems with Kp < 1 the essen- 
tial properties from the perturbative bosonic renormalization-group calculation and 
the boundary conformal field-theory analysis can be summarized as follows. The 
backscattering amplitude generated by a weak impurity is a relevant perturbation 
which grows as A^^p~^^^^, for a decreasing energy scale A, where z is the number 
of spin components. This behavior can be traced back to the power-law singularity 
of the 2/ci? density response function in a Luttinger liquid. On the other hand, the 
tunneling amplitude through a weak link between two otherwise separate wires is 
irrelevant and scales to zero as A"^, with the boundary exponent 

«B = -{K;' - 1) (2.3) 

depending only on the interaction strength and band filling, but not on the impurity 
parameters. At low energy scales any impurity thus effectively "cuts" the system into 
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2 Impurities in Luttinger liquids 



two parts with open boundary conditions at the end points, and physical observables 
are controlled by the open chain fixed point. 

In particular, the local density of states near the impurity is suppressed as 



p{uj) ~ \uj 



Ob 



(2.4) 



for \uj\ 0. 

Long-range Friedel oscillations in the density pr ofile induced by boundar ies or 
impurities decay with a power law at long distances |Egger and Grabertlll995l | as 



nix) 



X 



(2.5) 



for spinless fermions, where x measures the distance from the impurity or boundary. 
For spin- 1 fermions Kp is replaced by {Kp + l)/2 in Eq. ()2.5p . 

The conductance through an infinite Luttinger liquid with a single impurity 
vanishes at low temperatures as 



G{T) 



(2.6) 



The conductance through a single impurity of variable strength can be collapsed onto 
a single curve by a one-parameter scaling ansatz. For resonant scattering at double 
barriers the distance between the two barriers and the detuning from resonance 
introduce additional scales and a more complex behavior is observed. The Lorentzian 
resonance line shape for noninteracting electrons is modified by the interaction, and 



for appropriate parameters t 
as a function o f temperature 



le conductance exhibits distinctive power-law scaling 



Kane anc 



Furusakil Il998l : iNazarov and Glazman 
1994 . 



2nm[ iPolvakov and Gornvi 



Fisheiill992bl: iFurusaki and Nagaosalll993at 



201121; 



Vufi fit al. 



Note that the above power laws are strictly valid only in the absence of two- 
particle backscattering. For spin-| fermions they are in general modified by log- 
arithmic corrections. The asymptotic behavior is universal in the sense that the 
exponents depend only on the properties of the bulk system, via Kp, while they 
do not depend on the impurity strength or shape, except in special cases such as 
resonant scattering at double barriers, which require fine-tuning of parameters. 

The asymptotic low-energy properties of Luttinger liquids with a single impurity 
are rather well understood. Universal power laws and scaling functions have been ob- 
tained by bosonization, conformal field theory, and exact solutions for the low-energy 
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2.3 Experimental realization 



asymptotics in special integrable cases Giamarchill2004l | . Numerical results from ex- 
act diagonalization and DMRG applied to the lattice model of spinless fermions with 
nearest-neighbor interaction co nfirmed the field theoretical scenario and the validity 



of the underlying assumptions Eggert and Affleck 



1992; 



Meden et al. 



1998j . These 



methods are however limited to lattice systems with about 1000 sites and do not al- 
low for a systematic analysis of the crossover between the weak and strong-impurity 
limit. Moreover, only a restricted set of observables can be evaluated with affordable 
computational effort. In this context the fRG provides a complementary technique 
for microscopic models of interacting fermions with impurities, which does not only 
capture correctly the universal low-energy asymptotics, but allows one to compute 
observables on all energy scales, providing thus also nonuniversal properties, and in 
particular an answer to the important question at what scale the ultimate asymp- 
totics sets in. That scale can indeed be surprisingly low, and the properties above 
it very different from the asymptotic behavior. 



2.3 Experimental realization 



The progress in the fabrication of artificial low-dimensional structures led to ad- 
vanced experimental verification of the theoretical predictions. We present a short 
list of the most promising systems and of the employed experimental techniques. 
For a detailed discussion and references to the most recent p ublications and review 
articles on the subject we refer to Ref. 



Schonhammei 



2004 1 



Strictly one- dimensional systems are a theoretical idealization, the coupling to 

as the coupling between several Luttinger liquids 



Giamarchi 



200J|. The coupling between the chains 



an experimental probe as wel 
is not completely understood 
in a strongly anisotropic three-dimensional compound leads to the development of 
long-range order at very low temperatures in the phase for which the algebraic de- 
cay of the corresponding correlation function of the single-chain Luttinger liquid is 
the slowest. In appropriate temperature and energy regimes Luttinger-liquid be- 
havior can be expected in several systems with a predominantly one-dimensional 
character, as highly anisotropic quasi one-dimensional conductors, organic conduc- 
tors like the Bechgaard salts, as well as inorganic materials, artificial quantum wires 
in semiconductor heterostructu res or on surface sub strates, carbon nanotubes, and 
fractional quantum Hall fluids Schonhammerl l2004j . In particular, single- wall car- 



17 



2 Impurities in Luttinger liquids 



bon nanotubes are expected to show Luttinger-liquid behavior with Kp ~ 0.2 — 0.3 



down to very low 


emperatures, desi 


Dite t 


he presence of two low-energy 
I- 


Ee-g-er and Goffolin 


1997; 


Kane et al.1 


Il997 



Experimental techniques used to verify Luttinger-liquid behavior involve mainly 
high resolu tion photoemission and transport measurements, in addition to optical 



properties Schonhammeiil2004l |. A careful analysis of experimental data indicating 



power-law behavior and signatures of spin-charge separation reveals partly incon- 
sistent interpretations. The discussion on the modification of the quantized value 
/h for noninteracting electrons in a single channel by the interaction to Kp (e^//i) 
indicates a sensitive dependence on the schema tization of the conta cts, a challeng- 
ing theoretical as well as experimental problem S chonhammerl l20 04| . Experimental 
results for cleaved-edge overgrowth quantum wires and carbon nanotubes indicate 
power laws of the conductance consistent with Luttinger-liquid behavior. In the last 
few years, ultracold gases in optical lattices have opened up an entirely new area of 
physics, where strong correlations can be studied with unprecedented flexibility and 
control of the parameters. Further work is necessary for clear experimental evidence 
of Luttinger-liquid behavior. 
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3 Functional RG technique: a short overview 



We review the functional renormalization- group approach for interacting Fermi sys- 
tems in the IPI version. Introducing an infrared cutoff K in the free propagator and 
differentiating the effective action with respect to A, an exact hierarchy of differ- 
ential flow equations for the IPI vertex functions is derived, describing the gradual 
evolution from the microscopic model Hamiltonian to the effective action as a func- 
tion of the continuously decreasing energy cutoff. We briefly discuss the relation to 
alternative formulations of the fRG approach. 



3.1 Introduction 



The renormalization-group is a powerful method in the study of low-dimensional 
Fermi systems, providing in particular a systematic and unbiased method to study 
competing instabilities and entangled infrared singularities at weak coupling. Early 
renormalization-group approaches for one- dimensional systems, combined wi th exact 
solut i ons of fixed-poin t models, have been a major source of physical insight |s61vom 



1979; 



Giamarchi 



200J|. From the renormalization-group point of view, the existence 



of the Luttinger liquid requires t he cancellation of con tributions to the flow of the 
two-particle vertex to all orders Metzner et al.Ml998l |. This is a one-dimensional 
phenomenon, in higher dimensions the interactions in general diverge and the flow 
in the fermionic variables breaks down indicating a possible opening of a gap in the 

fermionic excitation spectrum. 

Wilson's renormalization-group approach Wilsonlll971 : Wilson and Kogutll974 | 



of successive integration of degrees of freedom with different energy scales deter- 
mines the evolution of the bare action of the system, given by the microscopic 
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3 Functional RG technique: a short overview 



Hamiltonian, to the final effective action, from whicli all physical quantities can 
be extracted. The hierarchy of coupled differential flow equations for the Green or 
vertex functions describing the full functional evolution of the effective action has 
been first implemented for boson i c field theories in the context of c ritical phenom- 



ena 



Wegner and Houghto: 



3 



1973; 



Polchinski 



1984; 



Wetterich 



1993j . The intuition 



of the relevance of fRG rnethods for interacting Fermi systems followed in th e 1990 s 



Benfatto and Gallavottilll990l : iFeldman 



and Trubowit? 



199"nl lShanka,Hn99ll 



Il994( . 



together with important rigorous work Salmhofen Il999l | . The infinite hierarchy of 
flow eq uations can be solv ed exactly only in special cases, for instance the Luttinger 
model Schiitz et al.ll2004| . Truncations however preserve the successive handling of 



energy scales and the consequent treatment of infrared singularities, characteristic 
of a renormalizat ion-group treatment. 

There are several variants of the fRG flow equations. The flow equations for 
the connected amp utated Green functions correspond to the Polchinski scheme. 



first derived in Ref. 



Po1chinskilh984 



j. The expansion of the con- 



nected am putated Green functions in IPI vertex functio ns led to the respective flow 
equations [Wegner and Houghton 1973 : Weinberg 1976|. subsequently derived frora 



Wetterich 



1993; 



Morrii 



1994 : 



the Legendre transform of the g enerating functional 

Salmhofer and Honerkamp boOl j. The Wick-ordered scheme is obtained from the 



Polchinski scheme by expanding the generating fun ctional of the connected ampu- 
tated Gree n functions in Wick-ordered polynomials Wieczerkowskilll988 : 



Salmhofei 



19981 Il999| . Important applications of the fRG in condensed-ma tter physics include 



the two-dirn ensional Hubbard model usin g the Polchinski scheme 



IM 1^ 



the Wick-ordered scheme 



Honerkamp et al 



Halboth and Metznei 



Zanchi and Schulz 



and also the 

200ll |. In the context of classical disordered systems 



IPI scheme 

a fRG appro ach is necessary to overcome the problem of dimensional reduction 



Wiese 



2003[. One-dimensional impurity problems and Luttinger-liquid physics are 



most conveniently investigated in the IPI scheme, as self-energy contributions are 
included to all orders. 

In the following the hierarchy of differential flow equations for the IPI vertex 
functions is derived, which is obtained by differentiating the corresponding generat- 
ing functional w ith respect to an infrared cutoff introduced in the free propagator 



Sa1mhofeT]h998 |. 
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3.2 Generating functional 



3.2 Generating functional 

We consider a system of interacting fermions with single-particle propagator of the 
noninteracting system Gq. The properties of the system are determined by the 
action 

^[^,^] = (^,GoV)-^[^,^], (3.1) 

where ip and ip are Grassmann variables associated with creation and annihilation 
operators, and is an arbitrary many-body interaction. Here we introduced 

the short-hand notation {ip, Gq^iP) = YIik k' i^K[Go^]K,K'ipK', where K contains the 
Matsubara frequency in addition to the single-particle quantum numbers and YIk 
stands for summation over the discrete indices and integrals over the continuous 
ones. 

All connected Green functions are obtained from the generating functional 
Negele and Orlandl 1987 1 defined by 



2o J 

dfiQltP,^] e-^[^'^] e-(^'^)-(*''^) , (3.2) 

with Grassmann source terms rj and f]. The normalized Gaussian measure with 
covariance Q = Gq ^ 

djjqlij, ^] = ^ diPdiP e('^"''3^) (3.3) 

includes the exponential of the quadratic part of the action and the noninteracting 
partition function Zq, such that / dfiQ[ip,ip] = 1. The generating functional for the 
connected Green functions is related to the partition function of the physical system 
with action ()3.1|1 by 

g[ri^f]] = -\nZ[ri,f]] . (3.4) 

In the noninteracting case = 0, and the Gaussian integral 

dfiQ [ij, i;] e-('^''')-(^''^) = e-(*''^°'') (3.5) 

implies that g[ri,f]] = (f/, Gq^)- 
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The connected m-particle Green functions are given by the derivatives of the 
generating functional Q[ri, f]] with respect to the source terms at 77 = f/ = 



Qm Qm 



dr]K[ ■ ■ ■ dr]K'^ dr]Km ■ ■ ■ 



Ti=r]=0 

(3.6) 



where (. . .)c is the connected average of the product of Grassmann variables between 
the brackets. 

The connected amputated Green functions are generated by the effective inter- 
action V[x, x] defined by 



e 



vix.x] = / ^^^[^^ ^] ^-v[^+x,i'+x] _ (3.7) 



The substitution x = GqI] and x = GqV^ where Gq is the transposed propagator, 
relates V[x, x] to the functional G[r],f]] by 

V[x,x] = G[v,v]-{r],Gov) . (3.8) 

The functional derivatives of V[x, x] generate connected Green functions divided by 
Go(-^i) • • • Go{Km) Go{K[) . . . Go{K'^), that is, propagators amputated from exter- 
nal legs in the corresponding Feynman diagrams. The term (f/, GqTi) cancels the 
noninteracting part of Q[ri, fj] such that V[x, x] = for V[ip, ip] = 0. Hence, the non- 
interacting propagator is subtracted from the one-particle Green function generated 

by V[x,x]- 

The generating functional r[0, 0] for the iP/ vertex functions 7^ is derived from 
the Legendre transform of G[f].,f]] by 

r[0, 4>] + (0, Q0) = g[r], n] + (0, v) - in, <P) , (3.9) 

with 

= ^ = — + 

df] (90 

(3.0) 
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3.3 RG differential flow equation for F 



and 



Q)-\ (3.11, 



5r] 6r] \ 5(f) Sep 

For the special case without interaction G[ri,f]] = (f/, Gq?]) leads to F[0, 0] = 0. 

The choice of the appropriate generating functional for a convenient formulation 
of a renormalization-group approach depends on the physical problem under i nves- 



tigat ion. For a detailed description of the different schemes we refer to Ref. Enss 
boost : ^^^^ ^ill concentrate on the IPI version of the fRG. 



3.3 RG differential flow equation for F 

In this section we briefly review the general renormalization-group setup, intro- 
duced as a transformation that leaves the generating functional for the correlation 
functions invariant, and concentrate subsequently on the derivation of a continuous 
renormalization-group equation for the IPI f unctions, following the derivation in the 



context of interacting Fermi systems in Ref. Salmhofer and Honerkampll2001 |. 

The addition principle for Gaussian fields implies that for the decomposition 
Go = Gq + Gq the corresponding Gaussian measure factorizes as 

= |d/XQ<[vl/<,^r<] y'd^Q>[vl/>,^r^]e-^[*<+*>+*'^<+^>+*] , (3.12) 

with ^ = ^< + The generating functional W corresponds to W = V for the 
particular choice y = V. This leads to the semigroup law of the renormalization 
group 

g-w[#,$] _ J ^^^^ [^/^ ^,j^„w>[vi.'+$,^'+4l ^ (3^^3) 

where in VV> = yV{Q^,y) the fields with propagator Q-^ have been integrated 
out. The semigroup law implies that the system {Q, y) under analysis is exactly 
equivalent to the system {Q^ ,W{Q^ ,y)). In the present case is a covariance 
with infrared cutoff A, and has support only for fields with energies smaller than 
A. 
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Set up in this way, the renormahzation group is simply a symmetry of the 
generating functional yV{Q,y). In differential form W{Q,y) is independent of A, 
that is, 

l.W{Q,y)=0. (3.14) 

Inserting the right-hand side of Eq. (|3.13|) leads to the flow equation describing 
the gradual evolution from y to the effective functional W as a function of the 
continuously decreasing energy cutoff A. >V> is an infinite power series in the fields; 
the quadratic and quartic terms correspond to the self-energy and the effective 
interaction, higher order terms are howe ver always present and the conver gence of 



the infinite series is a nontrivial problem Salmhofer and Honerkampll2001 



In the following the differential equation for the generating functional r^[0,0] 
of the IPI functions, starting point for the hierarchy of differential flow equations, 
is derived. Introducing an infrared cutoff at an energy scale A > in the bare 
propagator leads to a A-dependent generating functional for the connected Green 
functions defined by 

e-^^['''^l = j dfiQA[tP, ij] e-^l^''^! e-('^'^)-(*''^) . (3.15) 

The original functional is recovered in the limit A ^ 0. Similarly the functional 
(j)] generating the IPI vertex functions is constructed with Gq replacing Go in 
Eqs. ()3.9I - IXTT|) . Differentiating the above Eq. ()3.15|) for G^[i],f]] with respect to A 
yields 

= -[Tr(G^g^) + AQA]e-^'['''^ (3.16) 

where the first term comes from the derivative of the normalization factor ()3.3j) , and 
Tr denotes the sum over all space-time indices. The functional Laplace operator Ag 
is defined as 



24 



3.4 Expansion in the fields and exact hierarchy of flow equations 



The flow of ^^[?7, 77] is then 



Tr (G^g^) - Tr 



dA 



Q 



5rj 5fj 



5r] 



5r] 



(3.18) 



Using the Legendre transform ()3.9|) the derivative of r^[ 



reads 



OA 



Tr 



Tr 



5r] 5rj 



5rj 5rj 



(3.19) 



leading to the exact renormalization-group equation 



dA ^ 



Tr 



Tr 



dA 



With the initial condition 

T^"[<P,^] = V[<P,^] 
Eq. (j3.2()|l determines the flow of uniquely for all A < Aq. 



(3.20) 



(3.21) 



3.4 Expansion in the fields and exact hierarchy of 
flow equations 

The renormalization-group equations for the IPI m-particle vertex functions 7^ are 
derived by expanding rA[</), 0] in Eq. ()3.20j) as a power series in the fields. The 
coefficients in the expansion of 0] determine 7^ by 



00 m 

(3.22) 



[m\, 

m=0 ^ ' Ki...K„, K[...K'^ j=l 
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Due to the antisymmetry properties of the Grassmann variables only antisymmetric 
vertex functions contribute. 

Similarly the second derivative on the right-hand side of Eq. ()H.2()|1 can be 
expanded. Separating the ^-independent part corresponding to the self-energy yields 



2rAr 



->A I -pAr 



<t>=(j)=0 



(3.23) 



where the remaining functional f "^[0, 0] is defined by Eq. (j3.22|) with indices m and 
j starting from 2. The second term on the right-hand side of the flow equation ()3.2fl|l 
then reads 



2rAr 



6'T 



A\-l 



-1 



(f^ + (GY')"' 

oo 



(3.24) 



with the full propagator defined via the Dyson equation {G^) ^ = (Gq ) ^ — S^. 
Introducing the single-scale propagator as 



S^ = G' 



d{G^ 



A\-l 



OA 



G' 



(3.25) 



the differential equation ()3.2()|1 for is 



— = Tr 
dA 



Tr 



dA 



1=0 



d 



G^-^[G^]-\G'' T^y 



dA 



oo 

^(-i)'Tr \s^r^ {G^r^y 



(3.26) 



The first term corresponds to a vacuum energy not entering the correlation functions, 
while the second one contains one- loop diagrams with (/ + 1) vertices connected 
by one single-scale propagator and / full propagators G^. The term linear in 
generates self-energy corrections. 

Inserting the components 7^ on the left-hand side and the components 7^ on 
the right-hand side of the flow equation ()3.26p for we obtain a system of equations 
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3.4 Expansion in the fields and exact hierarchy of flow equations 



for 7^. In a graphical representation the equations for m < 3 are 




(3.27) 





(3.28) 






(3.29) 

The initial conditions for the vertex functions at A = cxo are given by the bare 
interactions of the system. In particular, the flow of the two-particle vertex starts 
from the antisymmetrized bare two-particle interaction while m-particle vertices of 
higher order vanish at A = oo, in the absence of bare m-body interactions with 
m > 2. 

Note that the right-hand side of the equation for 7^, contains 7m+i- The infinite 
system of differential equations contains only one-loop terms in every equation, as 
the differential formulation of Eq. ()3.2m) contains only a single trace, and for the 
IPI scheme no tree terms appear. The infinite hierarchy produces the full Green 
functions, generating graphs with an arbitrary number of loops; truncations amount 
to a partial inclusion of higher order contibutions generated during the flow, where 
the internal lines contain only r nodes above the cutoff scale A. Consequences of 
symmetries are discussed in Ref. Salmhofer and Honerkampll2001 |. 
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3.5 Comparison to other RG schemes 



Infrared divergencies arising in the context of perturbative expansions or in prox- 
imity of phase transitions can alternatively be regularized by temperature, a weak 
coupling strength or a finite system size. In the fRG approach the cutoff scale is 
introduced only in the quadratic part of the bare action, and the regularization is im- 
plemented with r espect to energy scales. The temperature and interaction f lows are 
derived in Refs. 



Honerkamp and Salmhofei 



temp- 

2nn4 



Honerkamp et al.M2004| respec- 
tively, a pedagogic introduction is given in Ref. lEnsslbnn^ The renormalization- 
group equations describe the flow of the correlation functions as the cutoff scale 
is lowered. The choice of the basis set for the correlation functions determines a 
particular scheme. 

In addition to the IPI scheme described previously, the various generating func- 
tionals introduced in Sec. 13.21 correspond to different schemes. Starting point for 
the Polchinski scheme is the effective interaction, generating functional of the con- 
nected amputated Green functions. The flow equation for V^[x, x] is derived by 
replacing Q by in Eq. ()3.7|) and taking the derivative with respect to A. An 
expansion in powers of x ^i-nd x of the functional V'^[x, x] iii the renormalization- 
group equ ation leads to Pq l chinski's flow equa tions for amputated connected Green 



functions 



Polchinskil 11984 iKeHer 



with a similar structure as for the 



connected Green functions. The connected amputated Green functions are the ex- 
pansion coefficients of the generating functional V^[x, x] terms of monomials of 
the source fields x X- Alternatively, one can also expand V^[x, x] with re- 
spect to Wick- ordered polynomials, leading to the Wick-ordered G reen functions as 
expansion coefficients [wieczerkowski 1988 : Salmhofei 1998 . [l999l |. The flow equa- 
tions are characterized by a bilinear structure in the vertices on the right-hand side 
connected by bare A-dependent propagators. The Wick ordering also implies that 
except for the differentiated propagator the internal lines are supported below scale 
A instead of above it. Thus, for a momentum cutoff only momenta close to the 
Fermi surface contribute at low cutoff scale A. This iustifle s a parametrization of 
the c oupling functions by projecting onto the Fermi surface Halboth and Metznei 



|. Self-energy corrections are however most conveniently taken into account in 
the IPI formalism with full propagators on the internal lines. 

In an exact treatment all schemes are equivalent, differences arise with trun- 
cations of the inflnite hierarchy of flow equations. While the full hierarchy of flow 
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equations leads to the correct solution to all orders in perturbation theory indepen- 
dently of the scheme, in the computation of the lowest orders a particular scheme 
might be more suitable than others, depending on the considered physical problem 
and properties. An important point for the choice concerns the possibility of an 
efficient parametrization of the effective interactions by a manageable number of 
variables. 

Continuous symmetries in the bare action lead to conservation laws and Ward 
identities relating Green and response functions, as a consequence of the Noether 
theorem. These are generally not preserved for the truncated flow equations, in 
cont rast to the solution of the inflnite flow-equation hierarchy, as shown in detail in 
Ref. |EnsJ2nn5| . For a gauge-invariant construction however, as for the temperature- 
flow scheme, the Ward identities between Green and response functions are satisfled 
exactly despite truncations. The related p roperty of self-consistency is satisfled by 
construction in conserving approximations Bavm and Kadanofllll961 |. but generally 
violated in truncated fRG flows. However, in the one- dimensional lattice models for 
Luttinger hquids, the truncated fRG is nevertheless surprisingly successful and self- 
consistency does not appear to play an important role. 
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4 Functional RG for Luttinger liquids 



We apply the fRG in the one-particle irreducible version to one-dimensional Fermi 
systems with impurities. The lowest order truncation of the fRG hierarchy of flow 
equations, where the two-particle vertex is approximated by the bare interaction, 
considered previously for spinless fermions, is extended including two-particle ver- 
tex renormalization, and generalized to spin-^ systems. For spinless fermions the 
quantitative accuracy of the results improves considerably, whereas for spin-^ sys- 
tems vertex renormalization is necessary to take into account that backscattering of 
particles with opposite spins at opposite Fermi points scales to zero in the low- energy 
limit. The underlying approximations are devised for weak interactions and arbitrary 
impurity strengths. Details on the computation of the relevant observables from the 
solution of the flow equations are presented. 

4.1 Microscopic models 

We consider various lattice fermion systems with spinless and spin-| fermions sup- 
plemented by different types of impurity potentials. The Hamiltonian has the form 

H = Ho + Hi + H,^p (4.1) 

where Hq is the kinetic energy, Hj a short-range interaction, and i^imp a static local 
or nonlocal impurity potential. 

We distinguish between spinless and spin-i fermions. 
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4.1.1 Spinless fermions 

For the spinless fermion model 

^ = -^E(4+iS + ^h+i) (4-2) 
i 

describes nearest-neighbor hopping processes with an amplitude t and 

Hi = UY^ rij Uj+i (4.3) 
j 

is a nearest-neighbor interaction of strength U, as shown in Fig. l4.1l We use standard 
second quantization notation, where cj and cj are creation and annihilation operators 
on site j respectively, and rij = is the local density operator. The impurity is 
represented by 

Hi^p = ^ Vj,, c], , (4.4) 
i.i' 

where Vjij is a static potential. For "site impurities" 

V,,, = V, (4.5) 
this potential is local. For the special case of a single site impurity 

Vj = V6,,, (4.6) 

the potential acts only on one site jo- We also consider "hopping impurities" de- 
scribed by the nonlocal potential 

^j'j = ^jf = ^i'j+i • (4-7) 

For a single hopping impurity 

hi-^i = it'-t)Sj,, (4.8) 

the hopping amplitude t is replaced by t' on the bond linking the sites jo and jo + 1- 
In the following we will set the bulk hopping amplitude t equal to one, that is, all 
energies are expressed in units of t. 

The clean spinless fermion model Hq + Hj is exactly soluble via the Bethe ansatz 



Yang and Yandll966l |. The system is a Luttinger liquid for all particle densities n 



31 



4 Functional RG for Luttinger liquids 



u t 




Figure 4.1: Spinless fermion model with nearest-neighbor hopping amplitude t, nearest- 
neighbor interaction U, and a local site potential V on site Jq. 



and any interaction strength, except at half filhng for \U\ > 2. For U > 2 a. charge 
density wave with wave vector vr forms; for U < —2 the system undergoes phase 
separation. The Luttinger-hquid parameter Kp, which determines all the critical 
exponents of t he liquid, can be computed exactly from the Bethe ansatz solution 
|Haldanjll98rJ |. At half filling Kp is related to U by the simple explicit formula 

K-^ = - arccos ( -- ) (4.9) 



P TT V 2 

for \U\<2. 



4.1.2 Spin-| fermions 

For spin- 1 fermions, the kinetic energy is given by 

^0 = + 4 w) ' (4-10) 



where Cj^ and Cj^ are creation and annihilation operators for fermions with spin 
projection a on site j. The interaction term of the extended Hubbard model contains 
a local interaction U and a nearest-neighbor interaction U' 

Hi = uY^ rij^ nji + U'Y^ Uj uj+i , (4.11) 
i i 

with rija = Cj^ Cj^ and nj = rij^ + riji, as shown in Fig. 14.21 For the pure Hubbard 
model only the local interaction U is finite, while U' = 0. The impurity term 

Himr> = J2J2^r,^W^j. (4.12) 
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U i t t 



Figure 4.2: Extended Hubbard model with nearest-neighbor hopping amplitude t, local in- 
teraction U and nearest-neighbor interaction U' . 



differs from the spinless case only by the spin sum. 

In the absenc e of impurities, the Hubbard model can be solved exactly using 
the Bethe-ansatz 



Lieb and Wu 



1968[, while the extended Hubbard model is not 



integrable. The Hubbard model is a Luttinger liquid for arbitrary repulsive in- 
teractions at all particle densities except ha lf filling, where the system becomes a 
Mott insulator |voitlll99,4 loiamarchil 120041 . The phase diagram of the extended 

ling it is a Luttinger liquid at 



Hubbard model is more complex. Away from half fi 
least for sufficiently weak repulsive interactions 



Voit 



1995|. For the Hubbard model 



the Lutti nger-liquid parameter Kr, can be computed exactly from the Beth e ansatz 



solution 



Frahm and KorepinlllQQOt iKawakami and YangHlQQO : 



Schuljh990| . 



4.2 Cutoff and flow equations 
4.2.1 Cutoff 

The cutoff introduced in Sec. 13.31 can be imposed in many different ways. The only 
requirement is that the infrared singularities must be regularized such that the flow 
equations allow for a regular perturbation expansion in powers of the renormalized 
two-particle vertex. Since translation invariance is spoiled by the impurity, a Mat- 
subara frequency cutoff is the most effici ent choice, while a mome ntum cutoff is less 



suitable. At T = the cutoff is sharp Andergassen et al.ll2004l |. the extension to 



T > will be addressed subsequently. The cutoff is imposed by excluding modes 
with frequencies below scale A from the functional integral representation of the 
system, or equivalently, by introducing a regularized bare propagator 

G^M = Q{\uj\ - A) GoM . (4.13) 
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Here Gq is the bare propagator of the pure system, involving neither interactions 
nor impurities. Instead of the sharp cutoff imposed by the step function G one 
may also choose a smooth cutoff function, but the sharp cutoff has the advantage 
that it reduces the number of integration variables on the right-hand side of the flow 
equations. Note that we will frequently write expressions which are well defined only 
if the sharp cutoff is viewed as a limit of increasingly sharp smooth cutoff functions. 
The suppression of frequencies below scale A affects all Green and vertex functions 
of the interacting system, which become thus functions of A. The original system is 
recovered in the limit A ^ 0. 



4.2.2 Truncation schemes 

The truncation of the fRG hierarchy of differential flow equations for the one-particle 
irreducible m-particle vertex functions and their parametrization with a manageable 
set of variables or functions leads to different approximation schemes. 

In the lowest order truncation of the fRG hierarchy (cf. Sec. 13.41) the renor 



malize d two-particle vertex is approximated by the bare interaction [Meden et al 
l2nn2atf h|. This truncation scheme, denoted by Scheme I, includes only the first 
equation in the hierarchy for the one-particle vertex function = — S^, where the 
self-energy T,^ is related to the interacting propagator by the usual Dyson equation 



G^= [(G^)-i -S^]"' . (4.14) 

Here and below G^, T,^ etc. are operators, which do not refer to any particular 
single-particle basis, unless we write matrix indices explicitly. The right-hand side 
of the flow equation for T,^ ()3.27|) involves the two-particle vertex and the single- 
scale propagator introduced in Eq. (j3.25j) which has support only on a single 
frequency scale Ic^l = A. The flow equation for the self-energy reads 

AsA(i', 1) = _ i ^ e*-^°^ S\2, 2') r'^(l', 2'; 1, 2) , (4.15) 

^ 2,2' 

where (3 is the inverse temperature. The numbers 1, 2, etc. are a shorthand for 
Matsubara frequencies and labels for single-particle states such as site and spin 
indices. Note that ui = uj[ and 1^2 = 1^2 due to Matsubara frequency conservation. 
The exponential factor in the above equation is irrelevant at any finite A, but is 
necessary to define the initial conditions of the flow at A = Aq — ^ 00. 
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For a sharp frequency cutoff the frequency sum on the right-hand side of the 
flow equation can be carried out analytically in the zero temperature limit, where 
the Matsubara sum becomes an integral. At this point one has to deal with products 
of delta functions — A) and expressions involving step functions 6(|ci;| — A). 
These at first sight ambiguous expressions are well defined and unique if the sharp 
cutoff is implemented as a limit of increasingly sharp broadened cutoff functions 
6e, with the broadening parameter e tending to zero. T he expressions can then 



be conveniently evaluated by using the following relation Morris! Il994| . valid for 
arbitrary continuous functions /: 

S,{x - A) f[Q,{x - A)] 6{x - A) [ f{t)dt, (4.16) 



where 6^ = Q'^. Note that the functional form of for finite e does not affect the 
result in the limit e — 0. 

In Scheme I in Eq. ()4.15j) is replaced by the antisymmetrized bare two-particle 
interaction r^,''2/.i2 — -^i',2';i,2, where the lower indices 1, 2, etc. label single-particle 
states (not frequencies). Since /i',2';i,2 is frequency independent, no frequency de- 
pendence of the self-energy is generated in the flow. Carrying out the frequency 
integration in the flow equation for the self-energy ()4.15|1 one obtains 

= E E ^'^"'^ GiA^^) /v,2';l,2 , (4.17) 
a;=±A 2,2' 

where 

G^{iuj) = [G^\iuj) - S'^]"' . (4.18) 

Note that has no jump at |ci;| = A, in contrast to G^. 

The flow is determined uniquely by the differential flow equation and the initial 
condition at A = oo. The self-energy at A = cxd is given by the bare impurity (site 
or hopping) potential V . In a numerical solution the flow starts at some large finite 
initial cutoff Aq. Here one has to take into account that, due to the slow decay of 
the right-hand side of the flow equation for S'^ at large A, the integration of the 
flow from A = oo to A = Aq yields a contribution which does not vanish in the 
limit Aq oo, but rather tends to a finite constant. Since G^^'i'''^) ~^ ^2,2' /ii^) for 
Ico"! = A ^ 00, this constant is determined as 

lim r dAY^Y. ^ = I E • (4-19) 
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Including the bare impurity potential Vi^y, the initial conditions for the self-energy 
at A = Ao ^ oo is 

= ^i,i' + ^5Z^i'.2;i.2- (4.20) 

2 

For the flow at A < Aq the factor e*"^"^ in Eq. ()4.17|1 can be discarded. 

A further development of the fRG approach for impurities in Luttinger liquids 
includes the two-particle vertex renormalization, denoted by Scheme II and used in 
the following if not specified otherwise. For spinless fermions this extension does not 
matter qualitatively, but the quantitative accuracy of the results improves consid- 
erably, in particular at intermediate interaction strengths. By contrast, for spin-| 
systems vertex renormalization is necessary to take into account that backscattering 
of particles with opposite spins at opposite Fermi points scales to zero in the low- 
energy limit. The right-hand side of the flow equation for the two-particle vertex 

()3.28j) involves itself, but also the three-particle vertex F^. Neglecting the 
contribution of the three-particle vertex to the flow of the two-particle vertex, the 
coupled system of flow equations for the two-particle vertex F'^ and the self-energy 
T,^ is closed. In terms of an expansion in the bare coupling function, this truncation 
is exact up to second order. However, the fRG provides more than just a second- 
order calculation: the evolution of the interaction and the self-energy is continually 
fed back into the fRG differential equation. This effectively sums up contributions 
from arbitrarily high orders and thus produces a scale-dependent resummation of 
perturbation theory. We note that it does not correspond to an expansion to a 
flxed loop order: the flow equations appear to be one loop, but they also take into 
account two-loop effects by iteration. The relevant question is whether higher or- 
ders signiflcantly change the flow, they certainly do so if the coupling functions get 
too large. The contribution of Fg to F'^ is small as long as F^ is sufficiently small, 
because Fg is initially (at Aq) zero and is generated only from terms of third order 
in F^. A comparison of the fRG results to exact DMRG results and exact scaling 
properties shows t hat the truncation er ror is often surprisingly small, even for rather 



large interactions [Meden et al.ll2002al lbl|. The explicit form of the truncated flow 



equation for the two-particle vertex reads 
(9 1 

- F^(l', 4'; 1, 3) F^(3', 2'; 4, 2) - (3 ^ 4, 3' ^ 4') 



^ F^(l',2';l,2)= i^^G^(3,3')5^(4,4')[F^(l',2';3,4)FV,4';l,2) 

3,3' 4,4' 
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+ r^(2', 4'; 1, 3) r^(3', 1'; 4, 2) + {3 ^ 4, 3' ^ 4') . 

(4.21) 

Diagrammatically, the individual contributions for the particle-particle and particle- 
hole channels written explicitly are 



2 2' 




1 1' 

1 1' 1 r 1 2' 1 2' 




2 2' 2 2' 2 1' 2 1' 

(4.22) 



Instead of solving the frequency integrated flow equation in full generality, we 
implement the following approximation: the frequency- dependent flow of the renor- 
malized two-particle vertex is replaced by its value at vanishing (external) fre- 
quencies, such that remains frequency independent. As a consequence, also the 
self-energy remains frequency independent. Since the bare interaction is frequency 
independent, neglecting the frequency dependence leads to errors only at second 
order (in the interaction strength) for the self-energy, and at third order for the 
vertex function at zero frequency. In addition to the quantitative errors we miss 
qualitative properties related to the frequency dependence of the self-energy, such 
as the suppression of the one-particle spectral weight in the bulk of a pure Luttinger 
liquid. On the other hand, a comparison with exact numerical results and asymp- 
totic analytical results shows that the impurity effects are not qualitatively affected 
by the frequency dependence of S, at least for weak interactions. 

The frequency-integrated flow equation for the two-particle vertex, evaluated at 
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vanishing external frequencies, has the form 



a;=±A 3,3' 4,4' 

+ ^3,3' (^^) G'4,4'(^^) {~^l',4';l,3 r3',2';4,2 + r2',4';l,3 r3',l';4,2) 

(4.23) 

with the initial condition Ty^^' i 2 — h',2';i,2- 

A crucial point is to devise an efficient parametrization of the vertex by a man- 
ageable number of variables. For a finite lattice system with L sites the fiow of the 
two-particle vertex T^, ^i-x 2 involves independent fiowing variables, if transla- 

tion invariance is assumed, and 0{L^) variables, if the infiuence of the impurity on 
the fiow of the two-particle vertex is taken into account. For a treatment of large 
systems it is therefore necessary to reduce the number of variables by a suitable 
approximate parametrization of the vertex. In the low-energy limit (small A) the 
fiow is dominated by a very small number of variables, t he others being irrelevant 
according to standard renormalization-group arguments jvoitl 1995 1. In particular. 



the frequency dependence of the vertex, discarded already above, is irrelevant for 
the fiow of at small A. For larger A one can use perturbation theory as a guide 
for a simple but efficient parametrization of F'^. 

We neglect the infiuence of the impurity on the fiow of the two-particle vertex, 
such that F^ remains translation invariant. While this is sufficient for capturing the 
effects of isolated impurities in otherwise pure systems, it is known that impurity 
contributions t o vertex renorma lization become important in macroscopically disor- 
dered systems |Giamarc'liill2nn4| . We also neglect the feedback of the bulk self-energy 
into the fiow of F"^, which yields higher order contributions in the renormalized in- 
teraction. The two-particle vertex is parametrized approximately by a renormalized 
static short-range interaction, which allows us to capture various features: the low- 
energy fiow of the vertex at fc^ in the pure system is obtained correctly to second 
order in the renormalized couplings; the nonuniversal contributions at finite energy 
scales are correct to second order in the bare interaction; the algorithm for the fiow 
of the self-energy remains as fast as in the absence of vertex renormalization, such 
that one can easily deal with up to 10'' lattice sites! 

For a more concrete treatment of the vertex renormalization, we now focus on 
a specific model. 
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4.2.3 Spinless fermions 

For spinless fermions the two-particle vertex and the self-energy are fully character- 
ized by either site or momentum variables. In the low-energy limit, the flow of the 
vertex is dominated by contributions with momenta close to the Fermi points, such 
that the right-hand side of the flow equation is determined by momentum compo- 



the vertex, there is only one such component which is nonzero: 



nents of the vertex F^, , with k^, fcg; ^ii ^2 — =t^F- E)ue to the antisymmetry of 

A?-|^ 



/ = ^tM,,-k, ■ (4.24) 

In the low-energy limit the momentum dependence of the vertex away from ^kp 
is irrelevant. There are therefore many possible choices for the functional form of 
F^, , , which all lead to the correct low-energy asymptotics. For a model with a 

^1 '^2'^1 '^2 

bare nearest-neighbor interaction U, a natural and efficient choice is to parametrize 
the flowing vertex simply by a renormalized nearest-neighbor interaction ^7"^, which 
leads to a real space vertex of the form 



F^' = U 

■II- 
Jl'J2'-?l'-?2 



with f/j^j2 = (<^ii,i2-i+^ii,j2+i) ■ This yields the following structure in momentum 
space: 

rf, , = 2f/^ [cos(A;; - k,) - cos{k'^ - k,)] 5^'"^ , , , (4.26) 

where the Kronecker 6 implements momentum conservation (modulo 27r) . The flow- 
ing coupling constant is linked to the value of the vertex at the Fermi points by 
the relation 

/ = 2U^ [1 - cos(2fci.)] . (4.27) 
The flow equation for (7^ becomes 

with the particle-particle and particle-hole contributions 
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—G^p-kp^^^) ^^p+kF^^^)^~kF,p+kF\kF,p-kF^V-kF,kF\p+kF-kF ' (4.29) 

where on the right-hand side of the flow equation is given by Eq. ()4.26p . Using 
Eq. (jOTj) to replace d^g^ by SaU^ on the left-hand side of Eq. (gSHI), one obtains 
a flow equation for U'^ of the simple form 

SkU^ = h{A) {U^f . (4.30) 

The function h{A) depends only on the cutoff A and the Fermi momentum kp. An 
explicit formula for h{A) can be obtained by carrying out the momentum integral 
in Eq. ()4.28|) using the residue theorem. For finite systems the momentum integral 
should be replaced by a discrete momentum sum; however, this leads to sizable 
corrections only for very small systems. Inserting the momentum structure of F^ 
()4.26|) into the flow equation ()4.28p and replacing by on the left-hand side 
yields 

^ = 9 • . E / ^ /(P' ^) ' 4.31 
aA 27rsm Ki;' Jo 

U) — dl A 

where 



2 sin^ /ci? sin^p (cos fci? — cosp)^ [cos(2fci?) — cosp] 
/(P'^) = 7— T^^^r—. 7^ 7-. 7^^. + 



{^u - - e-p) (^^ - i'p? {^u; - e^..,)(^^ - ^^.J ' 

(4.32) 

Here = —2 cos k — fiQ, with /ig = —2 cos kp, is the bare dispersion relation relative 
to the bare Fermi level. Since f{p,uj) can be written as a rational function of cosp 
and sinp, the p -integral can be carried out analytically using the substitution z = e*^ 
and the residue theorem. After a lengthy but straightforward calculation one obtains 
the following result for the coefficient h{A) in ()4.30p : 



h{A) = - Re 



i (/^o + ^A) 



2 ^"'^ ' ^ y (/io + zA)2 

3z/i^ - 10/igA - 12z/i2(A2 + 1) + 6AVo + 18A^o + GiA^ + ^A^ 
7r(2/io + ^A)(4 -ijI + A^- 2iAfio)^ 



(4.33) 
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The flow equation ()4.30|) can be integrated to 

where H{h.) is the primitive function of h{h) with H{A) 
h{A) one obtains 



(4.34) 

for A — > oo. Integrating 



+ 



-Re 

TT 



(4 - ^2)A2 _ 2^/io(2 - fil)A + 4-64 + 8 



2(4 - /i2)3/2 

^ sinh- ^ 



tanh 



2(4-/i2)^A2-22/ioA + 4- 
4 + /io + i/ioA 



v/(4 + /ig + i/ioA)2 + 4(A - 2i/ic 



where sinh~^ and tanh~^ denote the main branch of the inverse of the complex 
functions sinh and tanh respectively. 

At half filling, corresponding to kp = 7r/2, the function h{A) is particularly 
simple 

6 + A2 



MA) 



1 
2^ 



A 



(4 + A2)3/2 



such that U reduces to 



u 



A 



U/{2n) 



(4.36) 



(4.37) 



In Fig. 14.31 we show results for the renormalized nearest-neighbor interaction as 
obtained from the flow equation at various densities n, for a bare interaction U = 1. 
While the renormalization does not follow any simple rule at intermediate scales A, 
all curves saturate at a finite value U* in the limit A — » 0, c orresponding to a finite 



g*, as expected for a Luttinger-liquid fixed point |Voitlll99^ 

Parametrizing by a renormalized nearest-neighbor interaction has the enor- 
mous advantage that the self-energy, as determined by the flow equation ()4.17|) . is a 
tridiagonal matrix in real space, that is, only the matrix elements and ^fj±i are 
nonzero. Inserting from Eq. ()4.25|1 into (j4.17|l . one obtains the following simple 
coupled flow equations for the diagonal and off-diagonal matrix elements: 



OA 



A 



u)=±kr=±\ 
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Figure 4.3: Flow of the renormalized nearest- neighbor interaction for the spinless 
fermion model, for U = 1 and various densities n. 



Asa 



u 



2tx ^ 



lUO] 



(4.38) 



Note that the self-energy enters also the right-hand side of these equations, via 
= (Gq^ — SA)~i. Since and Gq^ are both tridiagonal in real space, the 
matrix inversion required to compute the diagonal and first off-diagonal elements 
of G^ from can be carried out very efficiently. An algorithm for the numerical 
solutio n of the flow equation for scaling linearlv with the system size is described 



in Ref. 



Andergassen et al .112004: 



Ensjl20o4 



Very large systems with up to 10'^ sites 



can be treated without extensive numerical effort. 



4.2.4 Spin-i fermions 



We now describe the parametrization of the s patial (or momentum) dep endences of 
the two-particle vertex for spin-| fermions Andergassen et al.ll2005bl | . employing 
a natural extension of our previous parametrization for the spinless case in Sec. 14.2.31 
We consider spin-rotation invariant lattice systems with local and nearest-neighbor 
interactions. This includes the extended Hubbard model. 

For a spin-rotation invariant system the spin structure of the two-particle vertex 
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can be decomposed in a singlet and a triplet part: 

T^ = TjS, , +TfT, , , (4.39) 

with 

S , , = — ( S ,5 1—5 ,5 , 

T , , = - ,5 , + 5 ,5 , V (4.40) 

Since the total vertex is antisymmetric in the incoming and outgoing particles, the 
singlet part has to be symmetric and the triplet part V^ antisymmetric. 

Proceeding in analogy to the case of spinless fermions in Sec. 14.2.31 we first list 
momentum components of the vertex with all momenta at ±A;^. For the triplet 
vertex the antisymmetry allows once again only one such component 

9t = ^t\kF-kp;kF-kF ■ (4-41) 

For the singlet vertex there are several distinct components at ifci?. Since we will 
neglect the influence of the impurity on the vertex renormalization, the renormal- 
ized vertex remains translation invariant. Hence the momentum components are 
restricted by momentum conservation: k'^ + = + ^2' modulo integer multiples 
of 27r. The remaining independent (not related by obvious symmetries) components 
are 

i/s2 s\kp ^—kp\kp ,—kp 

dti = '^s^\kp,kp;kp,kp (4.42) 

and in the case of half filling, for which kp = vr/2, also 

9% = r^7r/2,7r/2;-7r/2 -7r/2 ■ (4.43) 

The labels 2, 3, 4 are chosen in analogy to the conventional g-ologj notation for one- 



dimensional Fermi systems S61vomlll979l |. In order to parametrize the vertex in a 



uniform way in all cases, we will include the umklapp component g^^ not only at 
half filling, but at any density. The effect on the other components is negligible for 
the range of interactions and fillings considered. 

Extending our treatment of the spinless case in Sec. 14.2.31 we now parametrize 
the vertex by renormalized local and nearest-neighbor interactions in real space. 



43 



4 Functional RG for Luttinger liquids 



For the triplet part, there is no local component, and only one nearest-neighbor 
component compatible with the antisymmetry, namely 

Ul"" = , (4.44) 

which has the same form as the nearest-neighbor interaction in the spinless case. 
Note that Tf!^ does not depend on j and is equal to Fh For the 

symmetric singlet part, there is one local component 

^ = (4-45) 
and three different components involving nearest neighbors: 

TT/A _ -pA 

— ^ s\j,j+l;j,j+l 

pA _ -pA 

^ s\j+l,j + l;j,j 

w^i" = ri+i,,;,,, • (4.46) 

For the Hubbard model, the bare vertex is purely local and the initial condition for 
the vertex is given by = 2U, while all the other components vanish. For the 
extended Hubbard model, U',^" = U^^" = U' are nonzero. 

The triplet vertex is parametrized by only one renormalized real space coupling, 
which leads to a momentum representation of the form 

F^ , , , , , , = 2Ul^ [cos(fc; -K)- cos(k' - k,)] 5^'") , (4.47) 

where the Kronecker 6 implements momentum conservation (modulo 27r) . The flow- 
ing coupling Uj.^ is thus linked in a one-to-one correspondence to the Fermi momen- 
tum coupling by 

= 2U[^ [1 - cos{2kF)] (4.48) 

as in the spinless case in Sec. 14.2.31 In the singlet channel we have found four 
real space couplings, that is, one more than necessary to match the three singlet 
couplings in momentum space, gg2, 9^3, Q^a- We discard the interaction W^, because 
it does not appear in the bare Hubbard model, where it is generated only at third 
order in [/, while the pair hopping appears already in second-order perturbation 
theory. Fourier transforming the remaining interactions yields the singlet vertex in 
fc-space 
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-Ptcos{k, + k,)]5^^-l^^ (4.49) 



from which we obtain a hnear relation between the momentum space couphngs g^2' 
9s3y dsi renormahzed interaction parameters U^, U'^^, : 

g^, = + 2[/f [1 + cos(2A;^)] + 2P,^ 

gfs = - 4f/f - 2P.^ 

= Ut + ^U'J" + 2P^ cos{2kF) . (4.50) 

The determinant of this hnear system is positive for aU kp, except for kp = and 
TT. Hence the equations can be inverted for all densities except the trivial cases of 
an empty or completely filled band. 

We can now set up the flow equations for the four independent couplings U^^, 
U^, U'g^, and P^ which parametrize the vertex. Consider the case T = first. 
Inserting the spin structure ()4.39p into the general flow equation for the two-particle 
vertex ()4.23p . and using the momentum representation for a translation invariant 
vertex, the flow equation for the singlet and triplet vertices F^, for a = s,t, can be 
written as 

^ r"i = E E / 1 ( + ™ + ™' ) . 

a;=±A b,b'=s,f^ 

with the particle-particle and particle-hole contributions 

"^^) ^b\k[,k'2;p,ki+k2-p '^b'\p,ki+k2-p;ki,k2 



pp 


^PP ^0, 
— '^a,bb' '-^pl 


[iuj) 


^0 

'^ki+k2-p 


PH 


^PH ^0/ 
— ^a,bb' '-^pl 


[iu) 


'^p+ki~k[ 


PH' 


r<PH' ^0 
— ^a,bb' ^p' 







)'^b\k'2,p+ki-k'2;ki,p'^b'\p,k[;p+ki-k'2,k2 ' (4-52) 

The coefficients Ca,66' are obtained from the spin sums as 



ripp 


= 1 


r<PP 
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~ ^s,ts 
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^pp 
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^s,bb' — ^s,bb' ^t,bb' ~ ^t,bb' ■ [^.OO) 

Note that we have neglected the self-energy feedback in the flow of F^, such that 
only bare propagators Go enter. On the right-hand side of the flow equation we 
insert the parametrization ()4.47j) for and ()4.49|1 for F^. The flow of the triplet 
vertex F^ , , is evaluated only for (/c'^, /ci, ^2) = {kp, —kp, kp, —kp) as in 
Eq. ()4.4ip . which yields the flow of g^, while the flow of the singlet vertex F^ , , 

is computed for the three choices of {k[, ^2, /ci, ^2) which yield the flow of g^2y 9s3y 
(7^4 corresponding to Eqs. ()4.42|) and ()4.43|) . Using the linear equations ()4.48p and 
()4.50|1 to replace the couplings by the renormalized real space interactions on the 
left-hand side of the flow equations, we obtain a complete set of flow equations for 
the four renormalized interactions Uj.^, U^, f/^^, and of the form 

dj,U^ = J2 K'a" (A) Ut, U^. , (4.54) 



where a = 1, 2, 3, 4 labels the four different interactions. The functions h'^,^„{A) can 
be computed analytically by carrying out the momentum integrals in Eq. ()4.51|) via 
the residue theorem, details are reported in App. lA.ll The flow equations can then 
be solved numerically very easily. For finite systems the momentum integral should 
be replaced by a discrete momentum sum; however, this leads only to negligible 
corrections for the physical observables presented in Sec. 14.41 

After computing the flow of the real space interactions, one can also calculate 
the flow of the momentum space couplings g^ by using the linear relation between 
the two. In the low-energy limit (small A) one recovers the one-loop flow of the 
q-ology model , the general effective low-energy model for one-dimensional fermions 
Solvom I1974 for details see Sec. IA.21 In addition, our vertex renormalization 



captures also all nonuniversal second-order contributions to the vertex at ±kp from 
higher energy scales. 

In Fig. 14.41 we show results for the renormalized real space interactions together 
with the corresponding momentum space couplings, as obtained by integrating the 
flow equations for the Hubbard model at quarter filling and T = 0. Note that 
the couplings converge to finite fixed-point values in the limit A ^ 0, but the 
convergence is very slow, except for the momentum space couplings g^^ and (7^4. 
This can be traced back to the familiar behavior of the so-called backscattering 
coupling g^j_ = ^ {g^2 ~ fl'f )> that is, the amplitude for the exchange of two particles 
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with opposite spin at opposite Fermi points. Backscattering is known to vanish 
logarith mically in the low-energy limit for spin-rotation invariant spin-| Luttinger 
liquids We emphasize that this logarithmic behavior is not promoted 

to a power law by higher order terms beyond our approximation. By contrast, the 
linear combination of couplings which determines the Luttinger-liquid parameter Kp 
converges very quickly to a finite fixed-point value. 
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Figure 4.4: Vertex flow for the Hubbard model at quarter filling (n = l/2j and U = 1; 
upper panel.- flow of the renormalized real space interactions, lower panel; flow of the 
momentum space couplings. 
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Due to the above parametrization of the vertex by real space interactions which 
do not extend beyond nearest neighbors on the lattice, the self-energy generated by 
the flow equations is frequency independent and tridiagonal in real space. Inserting 
the spin and real space structure of into the general flow equation for the self- 
energy ()4.17|) . one obtains 



1 

Ait 



uj=±A 



r=±l 



47r 



AlUO] 



ui=±A 



lU] 



(4.55) 



where = (Gq ^ - S^)-^. 

Due to the slow decay of G at large frequencies, the integration of the flow 
equation for S from A = oo to A = Aq yields a contribution which remains finite 
even in the limit Aq oo, as described in Sec. 14.2.21 For the extended Hubbard 



model this contribution is given by S 



^1,1 



Ao 



U/2 + 2U' for J 



L — 1 and 



U/2 + U'. The numerical integration of the flow is started at a 



sufficiently large Aq with T,^° as initial condition. 



4.2.5 Extension to finite temperature 

At finite temperatures the Matsubara frequencies ciJ„ are discrete. The sum over 
Un of a function / can be written as an integral over a continuous variable u by 
introducing the distribution function P with a normalization f^^_^ |< duj P{uj) = 1 
for all n, 

E /K) = E / ^^^(^) = / du;P{u;)f{un,^) , (4.56) 

J |a;-a;n|< ttT J 

where a;„^^ denotes the discrete Matsubara frequency closest to u. Introducing a 
sharp frequency cutoff in the continuous variable u, the extension o f the flow equa- 



tions t o flnite temperatures is fairly simple, as pointed out by T. Enss |Andergassen et al, 
l2005bt . The general form of the flow equation for the generating functional for the 
IPI vertex functions T at T = 0, 

-^T^ = I dco5{\uj\-A)J^[ei\uj\^A),T\co)] , (4.57) 
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is modified to 



(4.58) 



at T > 0. Applying tlie lemma ()4.16|) the integral over uj can be carried out analyt- 
ically 



J^T^ = T j dujP{uj)6{\uj\-A) jytJ^[t,T^iur,,^)] 



(4.59) 



where F(-) = dt J^{t, ■). The contribution to the flow on the interval tUn, a — ttT < 
A < LJn,A + ttT is described by an autonomous differential equation, as the only 
explicit A dependence appears in P. As a consequence the result is independent 
of the particular choice of the distribution function P; for simplicity we choose the 
constant 



Piuj) 



2nT 



This leads to the final form of the flow equation 
d 



dA 



^7 



1 

2^ 



W=±IJ„, A 



[UJ 



(4.60) 



(4.61) 



for T. Hence, in the flow equations for the self-energy and the two-particle vertex, 
Eqs. ()4.17j) and ()4.2Hj) respectively, the expression u = ±A at T = is replaced by 
UJ = ±uJn, A at finite temperature, the functional dependence on oo remains the same. 
Note that for P{uj) = 6{lj — u;„_^) the flow equa tion cannot be s i mplified by 



()4.16|) and a smooth frequency cutoff has to be chosen [Enss et al 



2005; 



Enss^OOf 



4.3 Calculation of Kp 

The Luttinger-liquid parameter Kp, which determines the critical exponents of Lut- 
tinger liquids, can be computed from the fixed-point couplings as obtained from the 
fRG. A relation between the fixed-point couplings and Kp can be established via 
the exact solution of the fixed-point Hamiltonian of Luttinger liquids, the Luttinger 
model. A comparison of the fRG result for Kp with the exact Bethe-ansatz result 
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for the bulk model (without impurity) serves also as a check for the accuracy of our 
vertex renormalization. Since the above simplified flow equations yield not only the 
correct low-energy asymptotics to second order in the renormalized interaction, but 
contain also all nonuniversal second-order corrections at ±.kp from higher energy 
scales, the resulting Kp is obtained correctly to second order in the interaction. 



4.3.1 Spinless fermions 

For spinless fermions, Kp is determined by the Luttinger model parameters g and 
vp as 



Kp=, r -^/;^"^'^) (4.62) 

where g is the interaction between left and right movers and vp the effective Fermi 
velocity of the model, that is, the slope of the (linear) dispersion relation, with a 
possible shift due to interactio ns between particles moving in the same direction 
((74-coupling) already included IvoitlllQQ^ We therefore need to extract g and vp 
from the fRG flow in the limit A ^ 0. In order to obtain Kp correctly to order f/^, 
it is sufficient to obtain vp correctly to linear order in U . 

The Luttinger model interaction g and the fixed-point coupling g* = ^k^^kp-kp ~kp 
from the fRG are not simply identical, in contrast to what one might naively ex- 
pect. To find the true relation between g and g*, one has to take into account that 
the forward scattering limit of the dynamical two-particle vertex is generally not 
unique (in the absence of cutoffs), and depends on whether momentum or frequency 
transfers tend to zero first. This ambiguity is well-known in Fermi-liquid theory, 
where it lea ds to the distinction bet ween quasi-particle interactions and scattering 



amplitudes |Negele and Orlandlll987l | , but is equally present in Luttinger liquids, for 
the same reason in all cases: the ambiguity of the small momentum, small frequency 
limit of particle-hole propagators contributing to the vertex function. In the dynam- 
ical limit, where the momentum transfer q vanishes first, the singular particle-hole 
propagators do not contribute. In Fermi liquids this limit yields the quasi-particle 
interaction. In the opposite static limit the frequency transfer u vanishes first and 
particle-hole propagators yield a finite contribution. In the presence of an infrared 
cutoff A > the forward scattering limit of the vertex function is unique, since the 
ambiguity in the particle-hole propagator is due to the infrared pole of the single- 
particle propagator. Hence T^^^_k^.kp,-kp is well defined. However, T^^ _f^^.,^^ _i^^ 
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and also its limit for A — > depend on the choice for the cutoff function. For a 
momentum cutoff, which excludes states with excitation energies below A around 
the Fermi points, particle-hole excitations with small momentum transfers q are 
impossible. Hence particle-hole propagators with infinitesimal q do not contribute 
to the vertex at any A > 0, such that F^^ -kp-kp ~kp converges to the dynamical 
forward scattering limit, which is simply given by the bare coupling constant g 
in the Luttinger model. For a frequency cutoff the particle-hole propagators with 
vanishing momentum and frequency transfer yield a finite contribution at A > 0, 
which tends to the static limit for A — 0. This can be seen directly by integrating 
J2uj=±A I '^^^^ ^ from infinity to zero. Hence the vertex ^kp-kp;kp-kp 

obtained from our frequency cutoff fRG tends to the static forward scattering limit. 

For the Luttinger model the static forward scattering limit of the vertex can 
be obtained from the dynamical effective interaction between left and right movers 
D{q, iu), which is defined as the sum of particle- hole chains 

D{q,iy) = g + gli^_{q,iu) gli%{q,iu) g + ■ ■ ■ = q ^. ^ — — , (4.63) 



where 



Ul(q,iiy)=±- ^ (4.64) 

is the bare particle-hole bubble for right (-I-) and left (— ) movers. Note that only 
odd powers of g contribute to the effective interaction between left and right movers. 
This effective interaction appears naturally in the exact solution of the Luttinger 
model via Ward identities 



Dzvaloshinskii and Larkin 



1973 



Metzner and Di Castro 



1993J . For the static limit one obtains 

lim D(q, 0) = ^-^ —r (4.65) 

which we identify with our fixed-point coupling g* as obtained from the fRG with 
frequency cutoff. Inverting this relation between g and g* we obtain 



2'KVp 

3 = — — 
9 



^/invF^ + {g*y] . (4.66) 



For spinless fermions the difference between g and g* appears only at third order in 
the coupling, but for models with spin the distinction becomes important already 
at second order. 
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Figure 4.5: Luttinger-liquid parameter Kp as a function of U at various densities (as in 
Fig. \4-^ for the spinless fermion model; the inset shows the difference between the fRG 
result and the exact Bethe ansatz result for Kp. 

The Fermi velocity vp can be computed from the (frequency- independent) self- 
energy in momentum space as 

vf = v'p + dkJ^klk^ , (4.67) 

where Vp = dk€k\kp is the bare Fermi velocity. The self-energy is computed from the 
flow equation ()4.38p . which can be rewritten in momentum space as 

d . f dp 1 - cos(fc - p) 

where C,p = — /i- The chemical potential fi has to be fixed by the final condition 
^kp + ^kp = 0, where kp = vrn depends only on the density, not the interaction. 
From the tridiagonal structure of S in real space, but also from the above expression 
it follows that T,^ has the form = + cos k. The functional flow equation 
for Sj} yields a coupled set of ordinary flow equations for the coefficients and 
b^, with initial conditions = U and b^° = 0. The momentum integrals can be 
evaluated analytically via the residue theorem, such that the remaining set of two 
coupled differential equations (with as input) can be easily solved numerically. 
The result for vp is correct at least to first order in U, but not necessarily to second 
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order, since our simplified parametrization captures the two-particle vertex correctly 
to second order only at the Fermi points. 

Inserting g and Vp into the Luttinger model formula ()4.(i2|l we can now compute 
Kp as a function of U and density for the microscopic spinless fermion model. In 
Fig. 14.51 we show results for Kp{U) for various fixed densities. A compari son with 



exact results from the Bethe ansatz solution of the spinless fermion model |Haldane 
ll98rJ | in the inset shows that the fRG results are correct to second order in U and 
the vertex renormalization scheme described above is very accurate. 

4.3.2 Spin-| fermions 

The Luttinger-liquid parameter Kp for spin-^ fermions is given by 



The coupling constants gp2 and (7p4 parametrize forward scattering interactions in 
the charge channel (that is, spin symmetrized) between opposite and equal Fermi 
points respectively. They are related to the bare singlet and triplet vertices of the 
Luttinger model by 

fi'p2 ^ s\kp , — kp\kp ,—kp ~l~ '^^t\kp ,—kp\kp ,—kp ) 

gp/^ = -^ls\kp,kp;kF,kp ■ (4.70) 

These bare vertices are identical to the dynamical forward scattering limits of the full 
vertex F^. On the other hand, the vertex F^ obtained from the fRG with a frequency 
cutoff yields the static forward scattering limit for A ^ (cf. Sec. 14.3.11) . For the 
Luttinger model, the static forward scattering limit for the vertex can be computed 
from the effective interactions Dp2{q, ip) and Dp4^{q, iu), which are defi ned as the sum 
over all particle-hole chains with the bare interactions gp2 and gp4 SolvomI 1979 1. 



The summation becomes a simple geometric series if one introduces symmetric and 
antisymmetric combinations gp± = gp4±gp2 and Dp±{q, iv) = Dp^^q, iiy)±Dp2{q, iu). 
The static limit of the effective interaction Dp±{q,iu) yields the relation 

9;± = 1 ^-^f V (4-71) 

1 -^P±/(7rt'F) 
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between the Luttinger model couplings gp± and the fixed-point couplings 
from the fRG with frequency cutoff. Inverting 1)4.711) one obtains 



(4.72) 



9*p^K'^vf) 



(4.73) 



The Fermi velocity Vp can be computed from the self-energy for the translation- 
invariant pure system as in the spinless case, using the momentum representation 
of the fiow equations ()4.55p . 

The results for Kp from the above procedure are correct to second order in the 
bare interaction for the Hubbard model and also for the extended Hubbard model. 
While the fiowing couplings g^2 ^"^^ 9t converge only logarithmically to their fixed- 
point values for A ^ 0, the linear combination (7^2 + "^9^ which enters Kp converges 
much faster. 




n 



Figure 4.6: Luttinger-liquid parameter Kp for the Hubbard model as a function of electron 
density; results from the fRG are compared to exact results from the Bethe ansatz; the 
upper curves are for U = 1 and the lower ones for U = 2. 

In Fig. 14.61 we show results for Kp for the Hubbard model as obtained from the 
fRG and, for comparison, from the exact Bethe ansatz solution Frahm and Korepin 
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IQQOt iKawakami and Yandll990l : ISchulall990l |. Details on the solution of the cor- 
responding integral equations are reported in App. El The truncated fRG yields 
accurate results at weak coupling except for low densities and close to half filling. In 
the latter case this failure is expected since umklapp scattering interactions renor- 
malize toward strong coupling, even if the bare coupling is weak. At low densities 
already the bare dimensionless coupling U/vp is large for fixed finite U, simply be- 
cause vp is proportional to n for small n, such that neglected higher order terms 
become important. Note, for comparison, that for spinless fermions with a fixed 
nearest-neighbor interaction the bare dimensionless coupling at the Fermi level van- 
ishes in the low-density limit. 
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Figure 4.7: Luttinger-liquid parameter Kp for the quarter- filled extended Hubbard model 
as a function of U' for U = 0.5 and 1; results from the fRG are compared to DMRG data 
and to results from a one-loop g-ology calculation. 



For the extende d Hubbard model F ig. l4.7l shows a comparison of fRG results for 
Eiima et al.l 120051. The fRG results are exact to second order 



Kp to DMRG data 

in the interaction and are thus very accurate for weak U and U'. Results from a 
standard one-loop ^f-ology calculation as described in Sec. lA.2l deviate quite strongly 
already for U' > 0.5. In the g-ology approach interaction processes are classified 
into backward scattering {gi±), forward scattering involving electrons from opposite 
Fermi points {g2±), from the same Fermi points {g4±), and umklapp scattering {g3±)- 
All further momentum dependences of the vertex are discarded. This is justified by 
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the irrelevance of these momentum dependences in the low-energy limit, but leads to 
deviations from the exact flow at flnite scales, and therefore to less accurate results 
for the flxed-point couplings. 

The flow of gli_, i = 1, ...,4, is plotted in Fig. 14.81 in the upper panel for the 
quarter-flUed Hubbard model with bare interaction U = 1, and for the extended 
Hubbard model with U' = U / \pl in the lower. The fRG result is compared to the 
result from a one-loop (^-ology calculation. The backscattering coupling van- 
ishes logarithmically in both expected for the Luttinger-liquid flxed point 



Voit 



1995[. For the pure Hubbard model the good agreement with ^f-ology results 



stems from the purely local interaction in real space, since in that case pronounced 
momentum dependences of the vertex develop only in the low-energy regime where 
the (7-ology parametrization is a good approximation. By contrast, for the extended 
Hubbard model momentum dependences of the vertex which are not captured by 
the (7-ology classification (except for small A) are obviously more important. A gen- 
eralization of the (yf-ology parametrization of the vertex to higher dimensions, which 
amounts to neglecting the momentum dependence normal to t he Fermi surface, is fre- 



quent 



2000; 



V used in one-loop fRG calculations in two dimensions 



Zanchi and Schulz 12000: 



Halboth and Metzner 
Honerkamp et al.ll2001t iKampf and KataninI 2003| ^ 



The above comparison indicates that this parametrization works well for the pure 
Hubbard model, but could be improved for models with nonlocal interactions. The 
parametrization of the vertex by an effective short-range interaction used here could 
be easily extended to higher dimensions, where it will probably yield more accurate 
results, too. The relevance of an improved parametrization of the vertex beyond 
the conventional ^f-ology classification has also been demonstrated in a recent fRG 



analy sis of the phase diagram of the half-filled extended Hubbard model Tam et al, 

The inclusion of the momentum dependence due to the nearest-neighbor inter- 
action component in real space on the right-hand side of the fiow equation for the 
two-particle vertex modifies the fiow of the couplings at intermediate scales, before 
reaching the regime where a gf-ology description at weak coupling applies. A small 
repulsive initial backscattering amplitude may renormalize to an effective attractive 
one. For negative the renormalization group scales to strong coupling, indicat- 
ing an instability of the model towa rds a differ ent ground state characterized by a 
gap in the spin excitation spectrum |voitlh994 In Fig. 14.91 the phase boundary for 
the Luttinger liquid and spin gap is shown as a function of n and t/', as obtained 
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Figure 4.8: Flow of vertex on the Fermi points (in g-ology notation) at quarter filling 
and U = 1; upper panel; Hubbard model, lower panel; extended Hubbard model with 
U' = U /V^; the fRG flow is compared to the one-loop g-ology flow; note that in the upper 
panel fRG and g-ology results almost coincide. 



from the fRG together with the result from a one-loop g-ologj calculation. The fRG 
results confirm the spi n gap phase at lo w densities found with numerical Quantum 
Monte Carlo methods jclav et aPllQQd ]. where the spin gap develops with increas- 
ing U' from U' = 0. At low densities and close to half filling the truncated fRG 
results are not meaningful, since renormalization towards strong coupling occurs in 
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these limits and neglected higher order terms become important, see also Fig. 14.61 
Close to half filling the spin gap opens with increasing U' from [/' = f//2, as for the 
one-loop (yf-ology calculation. A full functional implementation of the momentum 
dependence of the two-particle vertex would allow a more detailed analysis of the 
phase diagram. 



1 - 



U = 1 



9i± > 



0.2 0.4 0.6 0.8 1 

n 

Figure 4.9: Phase boundary between the Luttinger liquid (gi± > 0) and spin gap phase 
(gi± < 0) for the extended Hubbard model as a function ofn and U' for U = 1; results from 
the fRG (solid lines) are compared to results from a one-loop g-ology calculation (dashed 
line). 



4.4 Observables 

In the next section we will present results for spectral properties of single-particle 
excitations near an impurity or boundary, the density profile and the linear con- 
ductance. Here we describe how the relevant observables are computed from the 
solution of the flow equations. 

4.4.1 Single-particle excitations 

Integrating the flow equation for the self-energy down to A = yields the physical 
(cutoff-free) self-energy S and the single-particle propagator G = (Gq ^ — From 



58 



4.4 Observables 



the latter spectral properties of single-particle excitations can be extracted. We focus 
on local spectral properties, which are described by the local spectral function 

Pj{uj) = ImGjjiu + iO+) , (4.74) 

where Gjj{uj + zO"*") is the local propagator, analytically continued to the real fre- 
quency axis from above. 

In our approximation the self-energy is frequency independent and can there- 
fore be viewed as an effective single-particle potential. The propagator G is thus the 
Green function of an effective single-particle Hamiltonian. In real space representa- 
tion this Hamiltonian is given by the tridiagonal matrix hcs = /lo + S, where the 
matrix elements of ho are the hopping amplitudes in Ho, Eq. ()4.2|1 . For a lattice with 
L sites this matrix has L (including possible multiplicities) eigenvalues eA and an 
orthonormal set of corresponding eigenvectors ipx- For the spectral function Pj{uj) 
one thus obtains a sum of 6 peaks 

P.M = 5^^A,<^(^-6), (4.75) 

A 

where = ex — /i, and the spectral weight wxj is the squared modulus of the ampli- 
tude of tpx on site j. For large L the level spacing between neighboring eigenvalues 
is usually of order L~^, except for one or a few levels outside the band edges which 
correspond to bound states. 

Due to even-odd effects etc. the spectral weight wxj generally varies quickly 
from one eigenvalue to the next one. A smooth function of u which suppresses 
these usually irrelevant finite-size details can be obtained by averaging over neigh- 
boring eigenvalues. In addition, dividing the spectral weight wxj by the level spacing 
between eigenvalues yields the local density of states, which we denote by Dj{uj). 



4.4.2 Density profile 

Boundaries or impurities induce a density profile with long-range Friedel oscillations, 
which are expected to decay with a power law at long distances |Egger and Grabert 
Il995| . The expectation value of the local density Uj could be computed from the local 
one-particle propagator Gjj, if G was known exactly. However, the approximate flow 
equations for S can be expected to describe the asymptotic behavior of G correctly 
only at long distances between creation and annihilation operator in time and/or 
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space, while in the local density operator time and space variables coincide. In the 
standard renormalization-grou p terminology rij is a composite operator, which has 
to be renormalized separately (Zinn-.TustinI 12002 1 . 

The flow equation for can be derived by computing the shift of the grand 
canonical potential generated by a small fleld (pj coupled to the local density. 
Alternatively to the numerical differentiation one may carry out the (pj derivative 
analytically in the flow equations, which yields a flow equation for the density in 
terms of the density response ve r tex. The g eneral structure at T = is described 



in Refs. 

the density is 

d 



Andergassen et al.ll2004l : lEnssll2005l |: the flnal form of the flow equation for 



1 

'2^ 



j=±A 



tr 



(4.76) 



with the density response vertex given by 



2,3 V^^) Rj;3,3 



' ^3',2'(^^) ri'^2';l,2 



(4.77) 



=±A 2,2' 3,3' 



Note that within the approximate treatment of the two-particle vertex described in 
Sec. 14. 2. 2^ the density-response vertex is frequency independent. 

For parametrized by local and nearest-neighbor interactions in real space, 
the matrix is tridiagonal, that is, only the components R^.^ 3 and Rf.3^3±i are 
nonzero. The initial condition for the density is = | for any fllling, due to 
the slow convergence of the flow equation ()4.76j) at large frequencies, which yields 
a flnite contribution to the integrated flow from A = 00 to Aq for arbitrarily large 
finite Aq, as in the case of the self-energy discussed in more detail in Sec. 14.2.21 The 
initial condition for the response vertex is Rflii = Sjidw. 

To avoid the interference of Friedel oscillations emerging from the impurity or 
one boundary with those coming from the (other) boundaries of our systems we 
suppress the influence of the latter by coupling the flnite chain to semi-infinite 
noninteracting leads, with a smooth decay of the interaction at the contacts, as 
described in detail in the next section. 



4.4.3 Conductance 

For the calculation of the conductance G a finite interacting chain (with sites 
1, . . . , L) is coupled to noninteracting leads at both ends, corresponding to an ex- 



60 



4.4 Observables 



perimental setup where the Luttinger-hquid wires are connected to (higher dimen- 
sional) Fermi-hquid leads. Using a proie ction rnethod the system with leads can 
be reduced to an effective L-site problem The leads are modeled by 

a one-dimensional tight-binding lattice with nearest-neighbor hopping amplitude t. 
The inffuence of the leads on the interacting chain can be taken into account by an 
additional dynamical boundary potential 



lead , 



iuJn + Ho 



4 



6 



iu}„) = ^— I i - w i - I [Oij + Olj) (4.78) 

The param- 



in the bare propagator Go of the interacting chain [Ens 
eter /ig is the chemical potential, which is related to the density n in the leads by 
fiQ = —2coskF. Uncontrolled conductance drops due to scattering at the contacts 
between leads and the interacting part of the chain can be avoided by switching off 
the interaction potential smoothly near the contacts. As long as the switching on of 
the interaction is smooth enough and the bulk part of the wire is large compared to 
the switching region, the results are independent of the microscopic details of this 



procedure. In addition, interaction-induc ed bulk shifts of the density have to be 



compensated by a suitable bulk potential |En; 




the current-current correlation function 



is comr 


)uted 


Mahan 


2nnn| 



20001 1 . Within our approximation 



scheme the self-energy h as no imaginary part, which implies that there are no vertex 
corrections jOgurilEoOl| . I n this approx imation the Kubo formula reduces to the 
Landauer-Biittiker formula [Pattal EqQbI . relating the transmission probability di- 
rectly to the linear conduct ance G(T,L) of noninteract ing fermions. For a detailed 
derivation we refer to Refs. 
conductance is given by 



Enss et al 



2005 



Ens£ 



2005l | ; the final expression for the 



G{T,L) = -z 



h 



\tie,T,L)\' fie) de 



(4.79) 



with |t(£,r,L)|2 = [4 - (/io + e)2]|Gi,L(£,r)|2, and / the Fermi function; z is the 
number of spin components. 
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In this section we present results for the local density of states near boundaries and 
impurities, the density profile, and the linear conductance for the spinless fermion 
model and the extended Hubbard model, as obtained from the solution of the fRG flow 
equations. A comparison with exact DMRG data is made for those observables and 
system sizes for which such data could be obtained. The asymptotic low-energy be- 
havior for weak and intermediate impurity strengths is approached only at rather low 
energy scales, accessible only for very large systems. For spin-^ systems two-particle 
backscattering leads to striking effects, which are not captured if the bulk system is 
approximated by its low- energy fixed point, the Luttinger model. In particular, the 
expected decrease of spectral weight near the impurity and of the conductance at low 
energy scales is often preceded by a pronounced increase, and the asymptotic power 
laws are modified by logarithmic corrections. 



5.1 Spinless fermions 

5.1.1 Effective impurity potential 

The typical shape of the self-energy representing the effective impurity potential can 
be seen in Fig. 15.11 where we plot the diagonal elements T,jj and the off-diagonal 
elements ^jj+i near a site impurity of strength V = 1.5 added to the spinless 
fermion model with interaction strength U = 1 at quarter filling. Recall that the self- 
energy is tridiagonal in real space and frequency independent within our treatment. 
The diagonal elements can be interpreted as a local effective potential, the off- 
diagonal elements as a nonlocal effective potential which renormalizes the hopping 
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amplitudes. At long distances from the impurity both T,jj and ^jj+i tend to a 
constant. The former describes just a bulk shift of the chemical potential, the 
latter a bulk renormalization of the hopping amplitude toward larger values. The 
oscillations around the bulk shifts are generated by the impurity. The wave number 
of the oscillations is 2kF = 7r/2, where kp is the Fermi wave vector of the bulk system 
at quarter filling. Fig. 15.21 shows S^^jo for a site impurity of strength V = 1.5 as a 
function of A for different system sizes L. For finite L the flow is effectively cut off 
on a scale ~ 1/L, a sequence of L provides an extrapolation to the thermodynamic 
limit. The renormalized potential at the impurity site remains finite in the limit 
L oo, while the Fourier transform S^^, for momenta with k — k' = ± 2kp 
diverges. Similarly for a hopping impurity the effective amplitude does not scale to 
zero in the limit L ^ oo, and the weak-link behavior is associated to the long-range 
oscillations in real space. The straight line in a log-log plot of the difference between 
the asymptotic value I]^^ j^{L 
power law dependence. 

0.5 



oo) and J2^^j^{L) shown in the inset indicates a 




Figure 5.1: Self-energy near a site impurity of strength V = 1.5 for the spinless fermion 
model at quarter filling and interaction strength U = 1; the impurity is situated at the 
center of a chain with L = 1025 sites. 

The amplitude of the oscillations generated by an impurity in S decays slower 
than the inverse distance from the impurity at intermediate length scales, but ap- 
proaches a decay proportional to — jo I for \j ~ jo\ ~^ oo. This can be seen most 
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T 1 — I—] 1 1 — 1 1 — 1 1 — 1 1 — 1 1 — r 




Figure 5.2: Sj^^ as a function of A for a site impurity of strength V = 1.5 for spinless 
fermions at half filling and U = 1; the impurity is situated at the center of a chain of 
length L; the inset shows the difference from the asymptotic value jo(°°) ~ ^io jo^^^ ^'^ 
a function of L. 



clearly by plotting an effective exponent /3j for the decay of the oscillations, defined 
as the negative logarithmic derivative of the oscillation amplitude with respect to 
the distance \j — jo I- In Fig- 15. HI we show the effective exponent resulting from the 
oscillations of Hjj as a function of the distance from a site impurity, for U = 1 and 
half filling. The impurity is situated at the center of a long chain with L = 2^^ + 1 
sites. To avoid interferences with oscillations from the boundaries we have attached 
semi-infinite noninteracting leads to the ends of the interacting chain, as described 
in Sec. 14. 4. HI Only for relatively large impurity strengths the asymptotic regime 
corresponding to /5j = 1 is reached before finite-size effects set in. For small V one 
can see that (3j increases from values below one, but the asymptotic long-distance 
behavior is cut off by the boundaries of the interacting region. For very small V (for 
example V = 0.01 in Fig. 15. Hj) we observe a plateau in Pj for intermediate distances 
from the impurity site. In this regime (3j is close to which can be un derstood by 
analytically solving the flow equations for small V Meden et al. 2002al lb|. 
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Figure 5.3: Effective exponent for the decay of oscillations of Sjj as a function of the 
distance from a site impurity of strengths V = 0.01, 0.1, 0.3, 1, 10 (from bottom to top), 
for the spinless fermion model at half filling and interaction strength U = 1; the impurity 
is situated at the center of a chain with L = 2^^ + 1 sites. 

5.1.2 Local density of states 

The long-range 2fci?-oscillations of the self-energy lead to a marked suppression of the 
spectral weight for single-particle excitations at the Fermi level, that is, at = 0. 
In Fig. 15.41 we show the local density of states Dj [u) on the site next to a site 
impurity of strength V = 1.5 for the spinless fermion model at half filling. The 
result for the interacting system at f/ = 1 is compared to the noninteracting case. 
Even-odd effects have been eliminated by averaging over neighboring eigenvalues 
(cf. Sec. 14.4.11) . 6 peaks outside the band edges corresponding to bound states are 
not plotted. The interaction leads to a global broadening of the band, which is due 
to an enhancement of the bulk hopping amplitude, and also to a strong suppression 
of Dj[uj) at low frequencies which is not present in the noninteracting system. For 
a finite system (here L = 1025) the spectral weight at the Fermi level remains 
finite, but tends to zero with increasing system size. In Fig. 15.51 we show results 
for the density of states choosing the same parameters as in Fig. EH but now for 
densities away from half filling: n = 1/4 and n = 3/4. In addition to the dip near 
u = a. second singularity appears at a finite frequency. This effect is due to the 
fact that a long-range potential with a wave number 2kF does not only strongly 
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scatter states with momenta near kp, but also those with momenta close to it — kp. 
Indeed the singularity is situated at = e^_A.> — where ek is the renormalized 
(bulk) dispersion. In the half-filled case only one singularity is seen simply because 
n — kp = kp for kp = 7r/2. 




Figure 5.4: Local density of states on the site next to a site impurity of strength V = 1.5 
for spinless fermions at half filling and U = 1; the impurity is situated at the center of a 
chain with 1025 sites; the noninteracting case U = is shown for comparison. 

Similar results are found for a hopping impurity, shown in Fig. 15.61 For an 
attractive interaction the density of states is strongly enhanced at the Fermi level. 
For comparison, typical results for the local density of states at a boundary are 
presented in Fig. 15.71 

The spectral weight at the Fermi level is expected to vanish asymptotically as 
a power law Ici^l"^, where 



as 



(5.1) 



is the boundary exponent describing the power-law suppression of th e density of 
states at the boundary of a semi- infinite chain with repulsive interactions Kane and Fisher 
That exponent depends only on the bulk parameters of the model, not on 



the impurity strength. For the spinless 
from the Bethe ansatz solution Haldane 



ermio n model it can be computed exactly 



1980 1 . 
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Figure 5.5: Local density of states on the site next to a site impurity as in Fig. \5.4\ (same 
parameters), but now for densities n = 1/4 and 3/4. 
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Figure 5.6: Local density of states on the site next to a hopping impurity t' = 0.25 for 
spinless fermions at density n = 0.6 and U = 1; the impurity is situated at the center of 
a chain with 1024 sites; the attractive case U = —1 is shown for comparison. 



In the above figures the envelope of the 6 peaks of weight w characterizing the 
spectral function of a finite system size introduced in Sec. 14.4.11 is shown. As a 
consequence the energy range over which a power-law suppression is observed is 
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Figure 5.7: Local density of states at the boundary, for the same parameters as in Fia. 15.61 

cut off by the finite size of tlie system at low energies. For a reliable analysis of 
the exponential behavior it is more convenient to consider the finite-size scaling 
of the spectral weight at the chemical potential. The large L dependence of the 
spectral weight is expected to exhibit a power law with the same exponent, and the 
scale where the power-law behavior sets in is given by L ~ uj/ttvf- In Fig. 15.81 the 
spectral weight at a boundary of a half-filled chain of length L = 10® and different 
interaction strengths U is shown as a function of u. The straight line in a log- 
log plot corresponds to a power law. Results for the same density and interaction 
parameters, but now as a function of system size L are presented in Fig. 15.91 Within 
the extension to finite temperatures (cf. Sec. I4.2.5j) . the same power-law behavior is 
found as a function of temperature, a detailed analysis in the context of transport 
phenomena follows below. 

We now analyze the asymptotic behavior of the spectral weight at the Fermi 
level by defining an effective exponent <y{L) as the negative logarithmic derivative 
of the spectral weight with respect to the system size, such that a{L) tends to a 
(positive) constant in case of a power law suppression. In Fig. 15.101 we show results 
for a{L) as obtained from the fRG for the spinless fermion model at quarter filling 
with up to about 10® sites, for a weak {U = 0.5) and an intermediate {U = 1.5) 
interaction parameter. The spectral weight has been computed either at a bound- 
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Figure 5.9: Spectral weight at the Fermi level as a function of system size L, for the same 
parameters as in Fig. \5.1A 



ary, or near a hopping impurity of strength t' = 0.5. Results obtained from the 
fRG without (upper panel) and with (lower panel) vertex renormalization, corre- 
sponding to Scheme I and Scheme II introduced in Sec. 14.2.21 are compared to exact 
numerical DMRG results (for up to 512 sites) and the exact boundary exponents 



69 



5 Solution of fRG equations and results 





Figure 5.10: Logarithmic derivative of the spectral weight at the Fermi level near a bound- 
ary (solid lines) or hopping impurity (dashed lines) as a function of system size L, for 
spinless fermions at quarter filling and interaction strength U = 0.5 (circles) or U = 1.5 
(squares); upper panel; without vertex renormalization, lower panel.' with vertex renorma- 
lization; the open symbols are fRG, the filled symbols DMRG results; the horizontal lines 
represent the exact boundary exponents for U = 0.5 and 1.5.; in the boundary case (solid 
lines) the spectral weight has been taken on the first site of a homogeneous chain, in the 
impurity case (dashed lines) on one of the two sites next to a hopping impurity t' = 0.5 in 
the center of the chain. 
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aB, plotted as horizontal lines. The fRG results follow a power law for large L, 
with the same asymptotic exponent for the boundary and impurity case, confirming 
thus the expected universality. However, the asymptotic regime is reached only for 
fairly large systems, even for the intermediate interaction strength U = 1.5. For the 
fRG Scheme I developed previously for impurities in spinless Luttinger liquids the 
effects of a single static impurity in a spinless Luttinger liquid are fully captured 
qualitatively, and in the weak-coupling limit also quantitatively. Originally devel- 
oped for the analysis of spectral densities of single-particle excitations, this scheme 
has been applied recently also to transport problems, such as persistent currents in 
mesoscopic r ings and the conductance of interacting wires con nected to noninter- 
acting leads [Meden et al.l 120031 : iMeden and Schollwockll2nn3al jbl]. The comparison 
with the exact DMRG results and exact exponents shows that the fRG Scheme II is 
also quantitatively rather accurate, and that the inclusion of vertex renormalization 
leads to a substantial improvement at intermediate coupling strength. 

A quantitative estimate of the accuracy of the exponents can be obtained from 
a comp arison to exact results from the Bethe ansatz solution of the spinless fermion 
model 



Haldane shown m Fig. 15.111 The results obtained without vertex 



renormalization (Scheme I) are represented with open symbols. The fRG results 
are correct to first order in the interaction U. In the approximate treatment of the 
two-particle vertex (cf. Sec. l4.2.2|l terms of order f/^ are only partially included, and 
an agreement to higher order can not be expected. Nevertheless the quantitative 
accuracy is improved considerably. The n dependence of the accuracy is related to 
the different importance of the vertex renormalization for different densities, as can 
be seen in Fig. 14.31 For half filling the vertex renormalization leads to a pronounced 
increase of the renormalized interaction, whereas for smaller fillings the effect is 
smaller, eventually leading to a decrease for n < 1/3. 

Results for the effective exponent a in the case of a site impurity are shown in 
Fig. 15. 12] at quarter filling and for an interaction strength U = 1. The comparison of 
the different curves obtained for different impurity strengths confirms once again the 
expected asymptotic universality, and also how the asymptotic regime shifts rapidly 
toward larger systems as the bare impurity strength decreases. The crossover scale 
depen ds on the bare impurity strength V as Kane and Fishej[l992al )~ 

E 



Fendlev et a1 



1 



V 



(5.2) 
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Figure 5.11: Difference between the fRG result and the exact Bethe ansatz result for the 
boundary exponent ub as a function of U , at densities n = 1/2 and 1/4 for the spinless 
fermion model as obtained from the power-law suppression of the spectral weight at the 
boundary; the open symbols are results without vertex renormalization, the filled symbols 
results with vertex renormalization. 

in agreement with our findings. The scale on which the impurity flows to strong 
couphng depends on the initial strength V and the "flow velocity" given by 1 — Kp. 
The effective flowing impurity strength can be estimated by 

Kfr ~ VL'-''^ , (5.3) 

where the crossover scale Lc in Eq. (|5.2|1 corresponds to Ves ^ irvp- Note that for 
weak impurities and intermediate system sizes the spectral weight follows a power 
law corresponding to the bulk behavior, and approaches the boundary exponent 
only at large distances. The bulk suppression of the spectral weight described by 
the anomalous dimension a = {Kp + — 2)/2 is not captured within the present 
scheme, since the self-energy is frequency independent. In the IPI version of the fRG 
the frequency dependence is generated by the two-particle vertex. The anomalous 
dimension could be included in an improved scheme by an iterative solution of 
the fRG flow inserting the two-particle vertex into the flow equation for the self- 
energy without neglecting its frequency dependence. This gives a two-loop diagram 
including the full two-particle vertex functions at scale A. The right-hand side of 
the differential equation for the self-energy is then nonlocal in A: the change of 
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the two-particle vertex and self-energy at scale A involves two-particle vertices and 
self-energies at scales A > A'. In presence of sufficiently strong effective impurity 
potentials however, the boundary behavior prevails over the bulk suppression of the 
spectral weight, as as ^ U and a ~ f/^. 
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Figure 5.12: Logarithmic derivative of the spectral weight at the Fermi level near a bound- 
ary (solid line) or site impurity (dashed lines) as a function of system size L, for spinless 
fermions at quarter filling and interaction strength U = 1; in the boundary case the spec- 
tral weight has been taken on the first site of a homogeneous chain, in the impurity case 
on the site next to a site impurity of strength V in the center of the chain; the horizontal 
line represents the exact boundary exponent for U = 1. 



5.1.3 Friedel oscillations 

We now discuss results for the density profile rtj. Boundaries and impurities induce 
Friedel oscillations of the local density with a wave vector 2k p- In a noninteract- 
ing system these oscillations decay proportionally to the inverse distance from the 
boundary or impurity. In an interacting Luttinger liquid the Friedel oscillations 
are expected to decay as |j — jol"'^'' at long distances |j — jo I- For a very weak 
impurity one expects a slower decay proportional to \j — jo\^~'^^'' at intermediate 
distances, and a crossover to the asymp totic power law with exponent Kp at very 



long distances |Egger and Grabertlll99q | . At intermediate distances the response of 
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the density to a weak impurity can be treated in linear response theory, such that 
the density modulation is determined by the density- density response function at 
2k F -I which leads to the power- law decay with exponent 2Kp — 1. 
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Figure 5.13: Amplitude (envelope) of oscillations of the density profile rij induced by 
a boundary as a function of the distance from the boundary, for spinless fermions with 
various interaction strengths U at half filling; the interacting chain with 2^^ + 1 sites is 
coupled to a semi-infinite noninteracting lead at the end opposite to the boundary; upper 
panel; log-log plot of the amplitude, lower panel: effective exponents for the decay, and the 
exact asymptotic exponents as horizontal lines. 

We analyze the long-distance behavior of the amplitudes more closely for the 
half-filled case, and compare to exact results for the asymptotic exponents. For 
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incommensurate filling factors the density profile looks more complicated, but at long 
distances from the boundary the oscillation amplitude has a well-defined envelope 
which exhibits a power law as a function of j. In Fig. 15.1 HI we show fRG results for 
the amplitude of density oscillations emerging from an open boundary, for a very 
long spinless fermion chain with 2^^ + 1 sites and various interaction strengths U. 
The end opposite to the open boundary is smoothly connected to a noninteracting 
lead. In a log-log plot (upper panel of Fig. l5.1Hj) the amplitude follows a straight 
line for almost all j, corresponding to a power-law dependence. Deviations from a 
perfect power law can be seen more neatly by plotting the effective exponent aj, 
defined as the negative logarithmic derivative of the amplitude with respect to j (see 
the lower panel of Fig. l5.1Hj) . The effective exponent is almost constant except at 
very short distances or when j approaches the opposite end of the interacting chain, 
which is not surprising. From a comparison with the exact exponent (horizontal 
lines in the figure) one can assess the quantitative accuracy of the fRG results. 



i ''//)■■' 




^ I I I I I I I I I I I I 

10^ 10^ 10=* 10-^ 10^ 

j - jo 



Figure 5.14: Effective exponent for the decay of density oscillations as a function of the 
distance from a site impurity of strengths V = 0.01, 0.1, 0.3, 1, 3, 10 (from bottom to top); 
the impurity is situated at the center of a spinless fermion chain with 2^^ + 1 sites and 
interaction strength U = 1 at half filling; the interacting chain is coupled to semi-infinite 
noninteracting leads at both ends. 

Effective exponents describing the decay of Friedel oscillations generated by site 
impurities of various strengths are shown in Fig. l5.14| for a half-filled spinless fermion 
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chain with 2^^ + l sites and interaction U = 1. Both ends of the interacting chains are 
coupled to noninteracting leads to suppress oscillations coming from the boundaries. 
For strong impurities the results are close to the boundary result (cf. Fig. 15. 11^ . as 
expected. For weaker impurities the oscillations decay more slowly, that is, with 
a smaller exponent, and approach the boundary behavior only asymptotically at 
large distances (beyond the range of our chain for < 1). For very weak impurities 
{V = 0.01 in Fig. I5.l4|l the oscillation amplitude follows a power law corresponding 
to the linear response behavior with exponent 2Kp — 1 at intermediate distances. 
Similar results are obtained for attractive interactions 



Andergassen et al 



2004 1 . 



A comparison with the exact expo nents from the Bethe ansatz solution of the 
spinless fermion model [HaldanJ ll98o| is shown in Fig. 15. 15^ where the difference 
AKp is reported as a function of the interaction U at densities n = 1/2 and n = 1/4. 
The open symbols represent the accuracy for the weak-impurity behavior in the 
linear response regime and the filled ones for the asymptotic exponent for strong 
impurities at long distances. Deviations from the exact result are quadratic in the 
interaction U. 




Figure 5.15: Difference between the fRG result and the exact Bethe ansatz result for 
the Luttinger-liquid parameter Kp as a function of U at densities n = 1/2 and 1/4 for 
the spinless fermion model, as obtained from the power-law decay of Friedel oscillations 
generated by a strong impurity at long distances (filled symbols), and by a weak impurity 
at intermediate distances (open symbols). 
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5.1.4 Scaling of the conductance 

In this section we study the transport through an interacting wire with a single 
impurity connected to semi-infinite noninteracting leads. The conductance exhibits 
the same asymptotic scaling behavior as a function of temperature T for an infinite 
wire as at T = as a function of L. For more than one impurity the T depen- 
dence is richer, showing nonmonotonic behavior and distinctive power-laws with 
differ ent universal exponents in va rious regimes; an extensive analysis is reported in 
Ref. Enss et al.ll2005t lEnsa l2005l |. Although this difference does not appear for a 



single impurity, we will mainly focus on the more physical temperature dependence 
of the conductance. 

Before analyzing the scaling behavior we present in Fig. a comparison 
of fRG results to numerical DMRG data for the conductance for a short wire at 
T = 0, determined from the persistent current observed in the presence of a m.a.g- 
netic fiux piercing a noninteracting ring in which th e interacting wire is embed- 
ded Meden et al. 2003 : Meden and SchoUwock 2003al )b|. The excellent agreement 



proves the reliability of the approximate fRG scheme for interactions in the range 
1/2 < Kp < 1. 
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Figure 5.16: Conductance as a function of the interaction U for a homogeneous spinless 
fermion chain at half filling, with L = 12 sites; the interaction is turned on sharply at the 
contacts. 



Fig. 15.171 shows typical fRG results for the T dependence of the conductance 
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through a single site impurity of strength V. The 1/T scahng observed for high 
temperatures (of the order of the bandwidth) results from the 1/T behavior of the 
derivative of the Fermi function in Eq. ()4.79|) . together with the weak temperature 
dependence of \t{e,T, L)\'^ at high T. For a strong impurity V = 10, G{T) follows 
a power law with exponent 2aB as indicated by the dashed line in Fig. 15.171 until 
saturation sets in for T ~ irvp/L. For an intermediate impurity the slope of the 
data tends towards the asymptotic exponent, but is still significantly away from it 
when finite-size saturation sets in. This slow change of the slope is a general feature 
of intermediate V. For a weak impurity G(T) approaches e'^/h. Similar behavior is 
found for the scaling of 1 — G(T) / (e^/ h) in the limit of a weak impurity predicted to 
follow T^'^^p"^^ which holds a s long as the correctio n to perfect conductance stays 
small 



Kane and Fisheiill992al ]d: iFendley et al.lll995 | 
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Figure 5.17: Temperature dependence of the conductance for a half-filled spinless fermion 
wire of length L = 10^, interaction U = 0.5 and a single site impurity of strengths V at 
the center of the wire; the dotted lines highlight power-law behavior. 



The above results are generic as long as the impurity is placed sufficiently away 
from the contact regions. The scale Sj^ = nvp/jo, where jo is the impurity position, 
sets a lower bound for the p ower- law scaling with the exponents discussed above 



id- 



Furusaki and Nagaosalll99(T| . For T ^ Sj^ a crossover to a power-law scaling with 



different exponents is found. For impurity positions in the contact region the ex- 
ponent as describes the tunneling between the noninteracting and the interacting 
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Luttinger liquids Enss et al.ll2005l |. Restrictions on the temperature range where 
universal scaling behavior might be detected arise from the bandwidth from above 
and the finite wire length from below. For an interacting wire of L lattice sites the 
energy scale 5l = nvp/L represents a lower bound for any temperature scaling. De- 
pending on the impurity and interaction strength an asymptotic low-energy regime 
might not be reachable in experiments on finite wires. For finite temperatures sys- 
tems of 10^ lattice sites are considered, comparable for typical lattice constants to 
quantum wires in the micrometer range accessible to transport experiments. 

Considering the conductance as a function of temperature and impurity strength, 
for a fixed Kp the renormalization-group flow from weak to strong impurity strength 
detern iines a scaling function G k ix) on which the data for different T and V col- 



lapse |Kane and Fisher ll 99230 
one-parameter scaling ansatz 



Moon et al .Ml 9931: iFendlev et al. 



199,4 



with 



x=[T/T^{U,n,V)] 



Ko-l 



Using a 



(5.4) 



the curves for G{T) and different V can be collapsed onto the i^'p-dependent scal- 
ing function Gxpix) for an appropriate nonuniversal scale TQ{U,n,V). It has the 
limiting behavior Gj^pix) ~ 1 — for x — > 0, and Gxpi^) ~ x~'^^^'' for x — * oo; 
for Kp = 1/2 and K„ = 1/3 the functional form of Gk^ was dete rmined explic- 
itly Kane and Fished Il992d : iMoon et al.l Il993l : iFendlev et al.l Il995l | . An example 
is shown in Fig. 15.181 for U = 0.5, the different colors stand for different impurity 
strengths V. As a consequence of the extended crossover region between weak and 
strong-impurity behavior, even for the fairly large system size of L = 10^ sites and 
the large range of temperatures we can treat, it is impossible to directly demonstrate 
the full crossover for a single set of parameters. A power-law behavior in both limits 
is found, with exponents which can be expressed consistently in terms of a single 
approximate Luttinger-liquid parameter Kp. Data for the same Kp but different 
interaction and filling parameters collapse on the results for half filling, since the 
scaling function depends on U and n only via the Luttinger-liquid parameter Kp. 
Considering different types of impurity potentials extending over more than one site 
or bond does not modify Gxp- 

A previous analysis of the zero temperature scaling behavior in the wire length 
L, replacing T in the above ansatz by ttvf/L, showed that one-parameter scaling is 
not affected by the presence of leads if the interaction is turned on very smoothly 
at the contacts and no one-particle scattering terms at the contacts are considered 
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Figure 5.18: One-parameter scaling plot of the conductance for the spinless fermion model 
at half filling with U = 0.5; the colors represent results obtained for different impurity 
strengths; the dashed lines indicate the asymptotic behavior for small and large x. 



Meden et al.ll2003l |. In addition, fRG data are found to collapse onto the K = 1/2 



Enss et al, 



local sine- Gordon scaling function known analytically |Meden et aLlboO^i 
l2nn5| . We finally remark that one-parameter scaling represents an excellent example 
for the power of the fRG technique capturing complex crossover phenomena at 
intermediate scales. 

There is an interesting correspondence between the transport through impuri- 
ties in a Luttinger liquid and the quantum Brownian motion in a cosine potential 
examined in Ref. 



Weiss 



The duality symmetry of the latter maps a weak 
impurity (small V) exactly to a strong impurity (large V) under the substitution 
Kp ^/ Kp- For Kp < 1 the system becomes localized, whereas the effective barrier 
height vanishes for Kp > 1. This symmetry can b e extended also to finite tem- 
peratures by a frequency- dependent transformation The equivalence 
between the mobility of the Brownian particle and the conductance through an im- 
purity in a Luttinger liquid relates the behavior for a strong impurity to the one for 
an appropriate weak impurity by 



atrong(T,i^p)/(eV/l) 



G'wcak(T,l/J^p)/(eV/i) . 



(5.5) 



In particular, the above relation holds for expansions around V <ti I and ^ 1. 
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The respective convergence radius defines the crossover scale between the two dual 
descriptions by T^^^°'^^{Kp) = TJ'^'^'^il/ Kp). For Kp > 1 high and low temperatures 
are exchanged. Similar scaling behavior is found for the nonl inear conductance, 



Weiss 



199 



where the bias assumes the role of the temperature 

As a consequence of this symmetry and the knowledge of the analytic form for 
a particular value of Kp ^ 1 the explic it scaling fu nction can be derived, as well 
as an expression for the crossover scale Weisj 19 gg|. The solutio n reproduces the 



result from the thermodynamic Bethe ansatz |Fendlev et al.lll995 |. 

5.2 Spin-| fermions 

5.2.1 Single-particle excitations 

For a; — * the spectral weights and the local density of states near a boundary 
or impurity are ultimately suppressed according to a power law with the boundary 
exponent 



with Kf, = 1 for spin-rotation invariant systems |Giamarchill2nn4| . However, due to 
the slow logarithmic decrease of the two-particle backscattering amplitude, the fixed 
point value of K„ is reached only logarithmically from above. Hence, we can expect 
that the asymptotic value of is usually reached only very slowly from below. 

The local density of states at the boundary of a quarter-filled Hubbard chain, 
computed by the fRG, is shown in Fig. l5.19l for various values of the local interaction 
U. Contrary to the expected asymptotic power-law suppression the spectral weight 
near the chemical potential is strongly enhanced. The predicted suppression occurs 
only at very small energies for sufficiently large systems. In the main panel of 
Fig. 15.191 the crossover to the asymptotic behavior cannot be observed, as the finite- 
size cutoff ~ irvp/L is too large. Results for a larger system with L = 10^ sites at 
[/ = 2 in the inset show the crossover to the asymptotic suppression, albeit only 
at very small energies. The dependence of the boundary spectral weight at the 
Fermi level on the system size L is plotted in Fig. 15.201 The L dependence of the 
spectral weight at zero energy is expected to display the same asymptotic power- 
law behavior for large L as the u dependence discussed above. Instead of decreasing 
with increasing L, the spectral weight increases even for rather large systems for 
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small and moderate values of U. For U > 2 the crossover to a suppression is visible 
in Fig. I5.2()[ For U = 0.5 only an increase is obtained up to the largest systems 
studied. The crossover depends sensitively on the interaction strength U, for small 
U it is exponentially large in nvp/U. 
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Figure 5.19: Local density of states at the boundary of a Hubbard chain of length L = 4096 
at quarter filling and various interaction strengths U ; the inset shows results for U = 2 
and L = 10^ at very low lo. 



The above behavior of the spectral weight and the density of states near a bound- 
ary of the Hubbard chain, that is, a pronounced increase preceding the asymptotic 
power-law su ppression, is captu r ed qualitatively already b y the Hartree-Fock ap- 
proximation Meden et al.l l2nnflt ISchonhammer et al.M200fll |. This is at first sight 



surprising, as the Hartree-Fock theory does not capture any Luttinger-liquid fea- 
tures in the bulk of a translational invariant system. In particular, a self-consistent 
Hartree-Fock calculation l eads to the unphy sical result of a charge- density- wave 
ground state for all f/ > |Cohen et al.lllQQsl . since a single impurity can not mod- 
ify bulk properties of the system. The initial increase of Dj{uj) near a boundary is 
actually obtained a lready within perturbation theory at first order in the interaction 



Meden et al 



2000| 



V{0) - zV{2ki 



■\n\uo/eF\ + 0{V'^ 



(5.7) 
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Figure 5.20: Spectral weight at the Fermi level at the boundary of a quarter-filled Hubbard 
chain as a function of system size L, for various different interaction strengths. 



where Dj{ujY is the noninteracting density of states, V{q) the Fourier transform 
of the real space interaction, and z the number of spin components. For spinless 
fermions {z = 1) with repulsive interactions the coefficient in front of the logarithm 
is always positive such that the first-order term leads to a suppression of Dj{uj). 
For the Hubbard model, one has z = 2 and V^(0) — 2V{2kp) = —U is negative for 
repulsive U. Hence, at least for weak U the density of states increases for decreasing 
uj until terms beyond first oder become important. For the extended Hubbard model, 
V^(0) — 2V{2kp) = 2U'[1 — 2cos(2fci?)] — U, which can be positive or negative for 
U,U' > 0, depending on the density and the relative strength of the two interaction 
parameters. At quarter filling V{0) — 2V{2kF) is negative and therefore leads to an 
enhanced density of states for U' < U/2. 

Using (yf-ology notation (cf. Sec. lA.2|) . one can write V{0)—2V{2kp) = g2±—'2gi±, 
which reveals that substantial two-particle backscattering {gi± > g2±/2) is necessary 
to obtain an enhancement of Dj{uj) for repulsive interactions. Backscattering van- 
ishes at the Luttinger-liquid fixed point, but only very slowly. In case of a negative 
V{0) — 2V{2kF) the crossover to a suppression of Djiu) is due to higher order terms, 
which are expected to become important when the first-order correction is of order 
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one, that is, for energies below the scale 

Uc = (-F exp , 

VV(0)-2V(2A;^) 

corresponding to a system size = ttvf/oJc- The scale uj^. is exponentially small 
for weak interactions. A more accurate analytical estimate of the crossover scale 
from enhancement to suppression h as been derived for the Hubbard model within 
Hartree-Fock approximation in Ref . Meden et al.ll2nnfll | . In a renormalization-group 
treatment uj^ is somewhat enhanced by the downward renormalization of backscat- 
tering. 

A comparison of fRG results with DMRG data for the spectral weight at the 
Fermi level is shown in Fig. 15. 2H for a boundary site in the upper panel, and near 
a hopping impurity of strength t' = 0.5 in the lower. The agreement improves at 
weaker coupling, as expected, and is generally better for the impurity case, compared 
to the boundary case. The deviations in the boundary case are probably due to 
our approximate translation-invariant parametrization of the two-particle vertex. 
Boundaries and to a minor extent impurities spoil the translation invariance of 
the two-particle vertex. Although the deviations from translation invariance of the 
vertex become irrelevant in the low-energy or long-distance limit, and therefore do 
not affect the asymptotic behavior, they are nevertheless present at intermediate 
scales. This feedback of impurities into the vertex increases of course with the 
impurity strength and is thus particularly important near a boundary. The scale for 
the crossover from enhancement to suppression of spectral weight discussed above 
depends sensitively on effective interactions at intermediate scales and can therefore 
be shifted considerably even by relatively small errors in that regime. 

With the additional nearest-neighbor interaction in the extended Hubbard model 
it is possible to tune parameters such that the two-particle backscattering amplitude 
becomes negligible. In that case the asymptotic power-law suppression of spectral 
weight should be free from logarithmic corrections and accessible already for smaller 
systems and at higher energy scales. The bare backscattering interaction in the ex- 
tended Hubbard model is given by gi± = U + 2U' cos{2kp) and therefore vanishes 
for U' = —U /[2 cos{2kp)], which is repulsive for U > if n > 1/2. In a one-loop 
calculation a slightly different value of U' has to be chosen to obtain a negligible 
renormalized g^j_ for small finite A, since the flow generates backscattering terms 
at intermediate scales even if the bare gi± vanishes. In Fig. 15.221 we show fRG and 
DMRG results 
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Figure 5.21: Spectral weight at the Fermi level near a boundary ('upper panelj and a 
hopping impurity t' = 0.5 ('lower panelj as a function of system size L for the Hubbard 
model at quarter filling and different interaction strengths U ; results from the fRG (open 
symbols) are compared to DMRG data (filled symbols). 

for the spectral weight of the extended Hubbard model at the Fermi level near 
a hopping impurity. In the upper panel a generic case with sizeable backscatter- 
ing is shown, while the parameters leading to the curves in the lower panel have 
been chosen such that the two-particle backscattering amplitude is negligible at 
low energy. Only in the latter case a pronounced suppression of spectral weight is 
reached already for intermediate system size, similar to the behavior obtained pre- 
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Figure 5.22: Spectral weight at the Fermi level near a hopping impurity t' = 0.5 as a func- 
tion of system size L for the extended Hubbard model with U' = U /V^, for various choices 
of U; upper panel; n = 1/2 (leading to sizeable backs cattering), lower panel.- n = 3/4 
(leading to small backscattering) ; results from the fRG (open symbols) are compared to 
DMRG data (filled symbols). 



viously for spinless fermi ons with nearest-neighbor interaction Andergassen et al 



20041 iMeden et al.ll2002al jbl|. This is also reflected in the energy dependence of the 
local density of states near the impurity. For parameters leading to negligible two- 
particle backscattering as in Fig. 15.231 the suppression of the density of states sets 
in already at relatively high energies and is not preceded by any interaction-induced 
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increase; the slight increase for small system sizes is a finite-size effect present also 
in the noninteracting case. Note also that the fRG results are much more accurate 
for small backscattering, as can be seen by comparing the agreement with DMRG 
data in the upper and lower panel of Fig. 15.221 especially for larger U. This indicates 
that the influence of the impurity on the vertex flow, which we have neglected, is 
more important in the presence of a sizable backscattering interaction. 



0.2 



0.1 - 



\'.' 



U = 
U = 0.5 
U = 1 
U = 2 



_L 



-0.5 




(jj 



0.5 



Figure 5.23: Local density of states near a hopping impurity t' = 0.5 in an extended 
Hubbard model with density n = 3/4 and interaction U' = U/V2 (leading to a small 
backscattering interaction) for various choices of U; the size of the chain is L = 4096. 



In the case of a negligible backscattering amplitude, the spectral weight at the 
Fermi level approaches a power law without logarithmic corrections for accessible 
system sizes if the impurity is sufficiently strong. The power law is seen most clearly 
by plotting the effective exponent a{L), that is, the negative logarithmic derivative 
of the spectral weight with respect to the system size. Fig. 15.241 shows a{L) on the 
site next to a site impurity of strength V for the extended Hubbard model with 
U = 1, U' = 0.65 and n = 3/4. The backscattering amplitude is very small for 
these parameters. The fRG results approach the expected universal ^-independent 
power law for large L, but only very slowly for small V. For weak bare impurity 
potential V, the crossover to a strong effective impur ity occurs only on a large 
length scale of order V^2/(Xp-i) Kane and Fishei 1992al lc|. For V = 0.1 this scale 



is obviously well above the largest system size reached in Fig. 15.241 The Hartree- 
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Fock approximation also yields power laws for large L, but the exponents depend on 
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Figure 5.24: Logarithmic derivative of the spectral weight at the Fermi level on the site 
next to a site impurity of strength V in the center of the chain as a function of system 
size L, for the extended Hubbard model with interaction parameters [7=1, U' = 0.65 and 
density n = 3/4; here the filled symbols are fRG, the open symbols Hartree-Fock results. 



The effective exponent obtained from the fRG calculation agrees with the exact 
boundary exponent to linear order in the bare interaction, but not to quadratic order. 
A quantitative estimate of the ac curacy is obtained from a comparison to exact 
DMRG results |Ejima et al.ll2005| shown in Fig. 15. 2 5^ for the extended Hubbard 
model at n = 3/4 and with U' = leading to a negligible backscattering 

amplitude. To improve this, the frequency dependence of the two-particle vertex, 
has to be taken into account. This is also necessary to describe inelastic processes 
and to capture the anomalous dimension of the bulk system (see also Sec. I5.1.2|l . 



5.2.2 Density profile 

Boundaries and impurities induce a density profile with long-range Friedel oscilla 
tions, which are expected to decay with a power law with expon ent (K^, + Kr,)l2 at 



long distances, where K^r = 1 for spin-rotation invariant systems Egger and Grabert 
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Figure 5.25: Boundary exponent as as a function of U at densities n = 3/4, with 
U' = U/\/2 (small backscattering interaction) as obtained from the power-law suppression 
of the spectral weight at the boundary; fRG results are compared to DMRG data. 



|l995[. For weak impurities linear response theory predicts a decay as \ j — 
at intermediate distances. 

As an additional benchmark for the fRG technique, we compare in Fig. 15.261 
fRG and DMRG results for the density profile Uj for a quarter-filled Hubbard chain 
with L = 128 lattice sites and open boundaries. Friedel oscillations emerge from 
both boundaries and interfere in the center of the chain. The fRG results have 
been shifted by a small constant amount to allow for a better comparison of the 
oscillations. Note that the mean value of nj in the tails of the oscillations deviates 
from the average density by a finite-size correction of order 1/L, which is related to 
the asymmetry of the oscillations near the boundaries. 

The long-distance behavior of the density oscillations as obtained within the 
fRG scheme has been analyzed in detail for spinless fermions in Sec. 15.1. HI For 
fermions with spin, asymptotic power laws can be identified only for special param- 
eters leading to negligible two-particle backscattering. In general, the asymptotic 
behavior of Friedel oscillations is realized only at very long distances, and the power 
laws are modified by logarithmic corrections. We finally remark the presence of 
a 4/cir-component of the Friedel oscillations for spin-i fermions, which decays as 
\j — iol"^^''- In present weak-coupling treatment this contribution is negligible. 
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Figure 5.26: Density profile nj for the Hubbard model with 128 sites and interaction 
strength U = 1 at quarter filling; fRG results are compared to DMRG data. 



since the 2kp component dominates for Kp > 1/3 Egger and Grabertlll995 |. 



5.2.3 Conductance 



Single impurity 

For a system of spinless fermions with a single impurity it was already shown that 
the conductance obtained from the truncated fRG obeys the expected power laws, in 



particular G{T) oc T^° ^ at low T, and one-parameter scaling behavior Meden et al 
2003 : Ieuss et ahlboO,^ . The corresponding scaling function agrees remarkably well 
with an exact result for Kp = 1/2, although the interaction required to obtain 
such a small Kp is quite strong. The more complex temperature dependence of 
the conductance in t he case of a doub l e barrier at or nea r a resonance is also fully 
captured by the fRG [euss et al .11200,4 Meden et a].ll2005 |. 

Fig. 15.271 shows typical fRG results for the temperature dependence of the 
conductance for the extended Hubbard model with a single strong site impurity 
{V = 10). Similar results were obtained for a hopping impurity. The considered 
size L = 10^ corresponds to interacting wires in the micrometer range, which is 
the typical size of quantum wires available for transport experiments. For U' = 
the conductance increases as a function of decreasing T down to the lowest tem- 
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peratures in the plot. For increasing nearest-neighbor interactions U' a suppression 
of G{T) at low T becomes visible, but in all the data obtained at quarter-filling 
the suppression is much less pronounced than what one expects from the asymp- 
totic power law with exponent 2a b- By contrast, the suppression is much stronger 
and follows the expected power law more closely if parameters are chosen such that 
two-particle backscattering becomes negligible at low T, as can be seen from the 
conductance curve for n = 3/4 and U' = 0.65 in Fig. 15.271 The value of Kp for these 
parameters almost coincides with the one for another parameter set in the plot, 
n = 1/2 and U' = 0.75, but the behavior of G{T) is completely different. Note that 
at T ~ TTVp I L finite-size effects set in, as can be seen at the low T end of some of the 
curves in the figure. An enhancement of the conductance due to backscattering has 
been found already e arlier in a renormalization-group stud y of impurity scattering 

JL 



in the (7-ology model 



Matveev et al. 



1993 



Yue et all 
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Figure 5.27: Temperature dependence of the conductance for the extended Hubbard model 
with L = 10^ sites and a single site impurity of strength V = 10, for a Hubbard interaction 
U = 1 and various choices of U' ; the density is n = 1/2, except for the lowest curve, 
which has been obtained for n = 3/4 and U' = 0.65 (leading to a very small backscattering 
interaction); the dashed line is a power-law fit for the latter parameter set. 

Results for the conductance of the extended Hubbard model with a hopping 
impurity with various amplitudes t' are shown in Fig. 15.281 The bulk parameters 
have been chosen such that the two-particle backscattering is practically zero at low 
T. From the plot of the logarithmic derivative of G{T) in the upper panel one can 
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Figure 5.28: Logarithmic temperature derivative of the conductance (upper panel) and 
of its deviation from the unitarity limit (lower panel) for the extended Hubbard model 
with L = 10^ sites and various hopping impurities. The density is n = 3/4, interaction 
parameters are U = 1 and U' = 0.65. The dashed horizontal lines highlight power-law 
behavior. 



see that for a strong impurity (small t') the conductance follows a well-defined power 
law G{T) oc T^°-s over a large temperature range. For intermediate t' the curves 
approach the asymptotic exponent at low T from below, but do not reach it before 
finite-size effects lead to a saturation of G{T) for T < irvp/L. For the weakest 
impurity in the plot, t' = 0.95, the conductance remains very close to the unitarity 
limit. However, the plot of the logarithmic derivative of 1 — G/{2e^ /h) in the lower 
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panel of Fig. 15.281 shows that 1 — G/{2e^ /h) increases as T^p~^ for decreasing T, as 
expected for a weak impurity in the perturbative regime Kane and Fisherlll992al ]c | . 
The effective exponents indicated by the two horizontal line s in the figure dey iate 
from the exact values (determined from the DMRG result Eiima et al.ll2005l | for 
Kp) by about 20% in the case of 2aB and only by 5% for Kp — 1 (cf. Sec. I5.1.4|) . 

Depending on the bare impurity and interaction parameters, nonuniversal be- 
havior dominates at intermediate energy and length scales. Moreover, in the pres- 
ence of backscattering the asymptotic power laws are modified by logarithmic cor- 
rections. 



Double impurity 

We fin ally present results for the conductance of a wire with a double-barrier im- 
purity Andergassen et al.ll2nn5aj . The setup modeling a quantum dot is shown in 
Fig. 15.291 Applying a gate voltage Vg on the dot sites j G [juir] by 



gate 



^9 '^i,- 



(5.9) 



the conductance can be tuned to resonance. Lo = jr — Jz + 1 is sufficiently far away 
from the contacts at sites 1 and L the position of the dot does not play a role. 
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Figure 5.29: Quantum dot schematization. 



Earlier studies of tunneling through a quantum dot embedded in a spinless Lut- 
tinger liquid showed that at T = and for finite L the resonances in the linear 
conductance G{Vg) are characterized by an almost Lorentzian shape with unitary 
height for symmetric barriers and a width w vanishing as a power law with a Kp- 



dependent exponent in the 
disappear for increasing L 



imit L 



Kane and Fishei 



oo. For asyrn metric barriers the resonances 
1992al l(j. At T > the peak value 
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of the conductance shows distinctive powe r-law behavior as a function of temper- 



ature 



Enss et ahl 120051 : iMeden et aLll2005l |. Including the spin degree of freedom 
the physics becomes more complex due to the appearance of the Kondo effect. For 
noninter acting leads (L = 1) the Ko ndo physics was investigated theoretically fo r 
the single-impurity Anderson model 



Glazman and Raikh 



198 



Ng and Lee 



Il988 |. 



At low temperatures and for sufficiently high tunnel barriers the Kondo effect leads 
to a broad plateau-like line shape of the resonance replacing the Lorentzian. On 
resonance the number of electrons on the d ot is odd imply ing a local spin-i degree 



of freedom responsible for the Kondo effect |Hewsonlll993l |. For the single- impurity 
Anderson model the conductance is p roportional to the one-par ticle spectral weight 
of the dot at the chemical potential Meir and Wingreenlll992| . Varying Vg within 
an energy range of order U the Kondo resonance of the spectral function i s pinned 
at Un at height 2e^ /h expl aining the broad plateau-like resonance in G{Vg) 



1993; 



Gerland et al. 



Hewson 



200 



^. The problem of a single s pin-^ coupled to a Luttinger 



liquid was investigated generalizing the Kondo model |Furusakill2005| . The fRG ap- 
proach allows for a direct computation of the electron transport through a quantum 
dot embedded in a Luttinger liquid in the presence of the Kondo effect. Here we 
address the question of the resonance line shape and the power-law scaling of G{T) 
resulting from the competition between the two correlation effects. 

For this purpose we first consider the situation L = L£) = la,tT = correspond- 
ing to the single-impurity Anderson model. Unless otherwise stated we consider 
symmetric dot-lead couplings. In Fig. 15.301 the conductance G as a function of gate 
voltage Vg for the single-impurity Anderson model is shown for different tunnel bar- 
riers t' in the upper panel, together with the oc cupation of the dot in the lower. For 
t' <^ U the resonance has a plateau-like shape (Gerland et al.ll2000l |. In this region 
the occupation is close to 1 while it sharply raises/drops to 2/0 to the left/right of 
the plateau. Also for asymmetric barriers we reproduce the exact resonance height 
4A2,Aij(AL + Ar)^ {2e^ /h), where A = Af, + Ap measures the hvbridizat ion of the 
dot and the (left and right) lead states Hewson 1993 : Gerland et allEooi^ . Here we 
are interested in the interplay of Kondo and Luttinger-liquid physics and thus focus 
on small t', that is, on tunnel barriers with small transmission. 

The pinning of the spectral function and the subsequent plateau-like resonance 
can be derived within the fRG Scheme I (cf. Sec. I4.2.2j) . For L = 1 the chemical 
potential only shifts the position of the resonance. For half filling G{Vg) is symmetric 
around and the fiow equation for the effective onsite energy V = Vg + ^fj^j^ on 
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Figure 5.30: Upper panel: conductance as a function of gate voltage for the Hubbard model 
at quarter filling, with U = 1, for L = = 1 and different t' ; lower panel: average 
number of electrons on the dot. 



the dot site jo reduces to 
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dA 



u 



vr 



Reg 



3D, 3D 



'^A) 



UV^/tt 



(A + A)2 + (yA)2 



(5.10) 



in the hmit of A ^ f/. Here A = 27rt'^p is the hybridization, and p den otes the 



Hewson 



spect ral weight at the end of the leads in the infinite band width limit 
Il993| . The initial condition is V^'^" = Vg. In this scheme the self-energy is frequency 
independent leading to a Lorentzian spectral function of width 2 A and height l/(7rA) 
centered around V = yA=o This implies that the spectral weight at po and thus 



G{Vg) is determined by V Meir and Wingreenlll992| | . The solution of the differential 



95 



5 Solution of fRG equations and results 



equation (j5.10|) at A = is obtained in implicit form 

vJi{v) - SJo{v) _ JojVg) 
vY,{v) - 6Yo{v) ~ Yoivg) ' 



(5.11) 



with V = Vn/U, Vg = Vgn/U, 6 = An/U, and Bessel functions J„, F„. For \Vg\ < Vc 
this equation has a solution with a small \V\, where Vc = Vcii/U is the first zero of 
Jo corresponding to Vc — 0.77U. For U ^ A the crossover to a solution with \V\ 
being of order U (for |V^| > V^) is fairly sharp. Expanding both sides of Eq. ()5.11|1 
for small \v\ and \vg\ gives 

^ = ^.^-p(^)- (5-12) 

The exponential pinning of the spectral weight at yUo = for small IV^^I and the 
sharp crossover to a of order U when \Vg\ > Vc leads to the observed resonance 
line shape. For U ^ A the width of the plateau is 2Vc — 1.5U, whic h is larger than 



thewidth U found with the numerical renormalization-group method Gerland et al, 



Our approximation furthermore slightly overestimates the sharpness of the 
box-shaped resonance. We expect the agreement is improved with the more accurate 
fRG Scheme II (cf. Sec. KT^ . 

In the complementary case of Luttinger-liquid leads and in the absence of Kondo 
effect, suppressed by turning off the interaction on the dot, similar results as for 
spinless fermions are obtained. Luttinger-liquid behavior leads to infinitely sharp 
resonances in the lir ait L — » cxj in str ong contrast to the broad resonances induced 



by the Kondo effect Enss et al.ll2005l |. To clearly observe Luttinger-liquid behavior 
for spin- 1 fermions at experimentally accessible scales one has to consider a situ- 
ation in which the backscattering process yielding logarithmic corrections to the 
power-laws is small by tuning the nearest-neighbor interaction U' , see Sec. 15. 2. II In 
Fig. I5.3ll the L dependence of GiVg) for a single-site dot computed for a small back- 
scattering amplitude and lead length L = 10^ typical for experiments is shown. At 
T = the width of the resonance in G{Vg) tends to zero for L —>■ oo. The extracted 
width w a.s a function of L reported in the lower panel with open symbols follows a 
power law lP^''~^^l'^ with an fRG approximation to the Luttinger-liquid parameter 
that is, correct to l eading order in the inte ra ction. Off resona nce G asymptotically 
vanishes as Kane and FisheJ 1992al )c: Ieuss et 311120051 . For Vg close to res- 



onance and 1 — G/{2e'^/h) -C 1, the deviation from the unitary limit increases as 
lJ--Kp characteristic of the scaling of a weak single impurity. Further increasing L 
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Figure 5.31: Upper panel: conductance as a function of gate voltage for the extended 
Hubbard model at n = 3/4, with U = 1, and U' = 0.65, for a noninteracting dot with 
Ljj = 1, t' = 0.1 and different L; lower panel: scaling of the resonance width for the 
same parameters as in the upper panel corresponding to a small backs cattering interaction 
(open symbols), and for n = 1/4, U = 1, U' = corresponding to a large backs cattering 
interaction (filled symbols). 



the behavior eventually crosses over to the off-resonance power-law suppression of 
G mentioned above. Due to an exponentially large crossover scale, even for the very 
large system sizes accessible with our method the complete crossover from one to the 
other p ower law can not be seen fo r a single fixed Vr , but follows from one-parameter 
scaling Kane and FisheJ 1992a c; Ieuss et all 120051 . For a sizeable backscattering 
amplitude the off-resonance conductance and thus the width first slightly increase 
for increasing L - becoming larger than for the noninteracting dot - followed by a 
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crossover to a decrease for exponentially large L, as indicated by the filled symbols in 
the lower panel. Due to the logarithmic vanishing of the backscattering process this 
behavior can in general not be observed on experimentally accessible scales, confirm- 
ing the important role of two-particle backscattering on intermediate length scales. 
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Figure 5.32: Upper panel: conductance G{Vg) as in Fia. \5.,yl\ (same parameters), but now 
with interaction on the dot; lower panel: scaling ofG/{2e^/h) atVg = Q outside the plateau 
(filled symbols), and of 1 — G/(2e^//i) on the plateau atVg = —0.685 (open symbols). 

We now analyze the linear conductance through a quantum dot in the presence 
of both Kondo effect as well as Luttinger-liquid leads. The upper panel of Fig. 15.321 
shows the L dependence of GiVg) for the same parameters as in Fig. 15. 3H but now 
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including the interaction on the dot. Note the different scales on the x-axis. For 
interactions large compared to the hybridization the broad plateau-like resonance 
induced by the Kondo effect is also present at least for finite Luttinger-liquid leads. 
The same holds for L^) > 1. The width of the plateaus is proportional to the 
ratio of the local component of the effective interaction at the end of the fRG flow 
and Lo- The differences between the curves for different L are barely visible, in 
particular the changes of the resonance width are marginal. For generic parameters 
with sizeable backscattering, the difference between curves computed for different L 
are even smaller. We note that the plateaus vanish if U, U', and n are chosen such 
that at the end of the fRG flow the local part of the interaction is small, and the 
resonance peaks are sharp. To analyze the L dependence at small backscattering 
in more detail in the lower panel of Fig. I5.H2I the scaling of G/{2e^/h) for a gate 
voltage outside the plateau (filled symbols) and of 1 — G/{2e^ /h) for a gate voltage 
on the resonance plateau (open symbols) are shown. Off resonance G follows a 
power-law with the exponent 2aB and G vanishes for L oo. Within every plateau 
there a value of Vg = VJ where G = 2e^//i independently of L. For Vg ^ VJ still 
within the plateau the deviation of G from the unitary limit scales as L^~^p , that 
is, with the weak single- impurity exponent. This shows that any deviation from VJ 
acts as an impurity. By analogy with the single-impurity behavior discussed in the 
previous sections we conclude that in the asymptotic low-energy limit the impurity 
will effectively grow and in the limit L ^ oo the plateaus will vanish. For infinitely 
long Luttinger-liquid leads the resonances are infinitely sharp even in the presence 
of Kondo physics. However, for tunnel barriers with small transmission the plateaus 
at finite L are well developed and the length scale on which the plateaus start to 
deteriorate is extremely large. For sizeable two-particle backscattering this scale is 
enhanced and the plateaus are more pronounced. 

Also for asymmetric barriers we find (almost) plateau-like resonances. To discuss 
this in more detail we focus on typical parameters with = 10^ and an asymmetry 
Ai/Ar ~ 2. Then the width is almost unaffected by the asymmetry. For the 
interaction and filling as in Fig. 15.321 the height within the plateaus varies by a 
few percent (with maxima at the left and right boundaries) and has average value 
~ 0.85 (2e^//i). With increasing L the variation of the conductance on the plateaus 
increases while the average value decreases. We expect that for L — > oo the resonance 
disappears. 

Concerning the power-law scaling of the conductance as a function of tempera- 
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5 Solution of fRG equations and results 




Figure 5.33: Temperature dependence of the conductance for the extended Hubbard model 
with L = 10^ sites, Hubbard interaction U = 1, for U' = 0.65 and n = 3/4 (leading to 
a small backs cattering interaction), and for U' = 0.75 and n = 1/4 (leading to a sizable 
backs cattering interaction); dot parameters: t' = 0.1, = 100; upper panel.' log-log plot 
of the conductance, lower panel.' effective exponents; the gate potential is chosen at the 
center of the resonance plateau closest to Vg = 0. 



ture, G{T), the resonant tunneling behavior in general extends over a range of gate 
voltages defined by the width of the resonance plateau, since for all experimentally 
accessible length scales and for typical asymmetries of the dot-lead hybridizations 
the plateau-like resonances characteristic for Kondo physics will also be present if 
the leads are Luttinger liquids. The vanishing of the resonance plateau in the limit 
of infinite system size is beyond the infrared cutoff scale set by the system size 
T ~ irvp/L for the appearance of power-law scaling with interaction-dependent ex- 
ponents characteristic for Luttinger-liquid beh avior. Hence the behavior of G(T) i s 
similar to the one found for spinless fermions Enss et al.ll2nn5l : iMeden et al.ll2nn5 |. 
In Fig. 15.331 the temperature dependence of the conductance for a large dot with 
Lfj = 100 and for Vg placed at the center of the resonance plateau closest to = 



100 



5.2 Spin-i fermions 



is shown. We distinguish the case of sizable and small backscattering processes 
9i±- -^o^ ^ small backscattering amplitude we identify temperature regimes in which 
G{T) follows distinctive power-law behavior with universal exponents. For T larger 
than the level spacing of the dot given by 7ivp/L£, a power law with the single- 
impurity exponent 2aB is found, arising from two independent single barriers acting 
as resistors in series. It follows a sequential tunneling regime characterized by the 
power-law exponent — 1, until saturation sets in at T ~ irvp/L. Small dots 
with L/) of order 1 exhibit only the latter. Similar results are obtained for other Vg 
within the resonance plateau, outside the plateaus the conductance follows the on- 
resonance behavior down to a scale set by the deviation from resonance. For smaller 
T we find a crossover to T^"^. By contrast, for a sizable backscattering amplitude 
no clear power law can be distinguished, as can be seen from the conductance curve 
for n = 1/4 and U' = 0.75 in Fig. 15.331 The value of Kp for these parameters almost 
coincides with the one for another parameter set in the plot, n = 3/4 and U' = 0.65, 
but the behavior of G{T) is completely different. For weak barriers analogous power- 
laws are detected in l — G/(2e'^/h), described by the single-impurity exponent Kp — 1 



at high temperatures and by Kp -|- 1 in the seque ntial tunne 



Enssetal 



ing re gime, similarly to 



2005 1 . 



the spinless case investigated in detailed in Ref. 

For systems with long-range interactions backscattering is strongly reduced 
compared to forward sca t tering. This seem s to be the case in carbon nanotubes 



Egger and Gogo1inlh997 



Kane et al.l Il997| . Hence, the conductance can be ex- 
pected to follow the asymptotic power law at accessible temperature scales for suf- 
: iciently strong impurities in these systems, as is indicated also by experiments 



Yao et al.lll999| . However, the effects due to two-particle backscattering should be 
observable in systems with a screened Coulomb interaction. Whether Luttinger- 
liquid behavior has convincingly been demonstrated experimentally is still a matter 
of debate. Nonetheless, the above scenario has to be taken into consideration in 
the discussion of the influence of boundaries or impurities in quasi one-dimensional 
conductors. 
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6 Conclusions and outlook 



The fRG provides a powerful tool to compute the intriguing properties of Lut- 
tinger liquids with static impurities. It captures the physics at all energy scales 
from the Fermi energy to the ultimate low-energy limit. The presented compu ta- 
tion scheme extends previous work for spinless fermions Meden et al.ll2nn2al ]bl] to 
spin- 1 fermions, including the vertex renormalization in addition to the renormaliza- 
tion of the effective impurity potential. The underlying approximations are devised 
for weak short-range interactions and arbitrary impurity potentials. The results 
agree remarkably well with exact asymptotic results up to intermediate interaction 
strength, and cover the universal low-energy asymptotics, as well as nonuniversal 
behavior and crossover phenomena at higher scales. 

Various observables have been computed for different fermion lattice models: 
spectral properties of single-particle excitations, the oscillations in the density pro- 
file induced by impurities or boundaries, and the temperature dependence of the 
linear conductance. The comparison to DMRG results, for those observables and 
system sizes for which such data could be obtained, yields a good agreement at weak 
coupling. For intermediate interaction strengths with sizable two-particle backscat- 
tering and strong impurities the deviations are significantly larger for spin-i than for 
spinless fermions. For the computation of the nonuniversal behavior at intermediate 
scales the neglected influence of impurities on vertex renormalization, in particular 
the interplay of impurities and the two-particle backscattering amplitude, is proba- 
bly important for fermions with spin. 

We confirm the universality of the open chain fixed point, but it turns out 
that very large systems are required to reach the fixed point for realistic choices of 
the impurity and interaction parameters. The spectral properties of single-particle 
excitations and the Friedel oscillations in the density profile induced by impurities 
or boundaries present the characteristic asymptotic power laws at low energy or 
large distance. For the linear conductance through a single impurity in Luttinger 
liquids connected to noninteracting leads the fRG captures the expected power-law 
scaling, as well as the complete crossover from the weak to the strong-impurity 
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limit determined by a one-parameter scaling function. For resonant tunneling in a 
Luttinger liquid with a double barrier enclosing a dot region, depending on the dot 
parameters several temperature regimes with distinctive power-law behavior of the 



resonance conductance as wel 
are identified 



Enss et al 



2005 



as regimes charact erized by non-universal behavior 



Meden et al 



200 



a 



Including the spin degree of freedom, two-particle backscattering of particles 
with opposite spin at opposite Fermi points leads to two important effects, not 
present in the case of spinless fermions. First, the expected decrease of spectral 
weight and of the conductance at low energy scales is often preceded by an increase, 
which can be particularly pronounced for the density of states near an impurity or 
boundary as a function of uj. For the density of states near a boundary this effect has 
been fo und already earlier within a Hartree-Fock and DMRG study of the Hubbard 
model 



Meden et al, 



2000; 



Schonhammer et al 



2000^, and for the conductance by a 



renor malization-group analysis of the gf-ology model 



Matveev et al. 



1993; 



Yue et al 



Second, the asymptotic low-energy power laws are usually modified by log- 
arithmic corrections. In the extended Hubbard model the backscattering can be 
eliminated for a special fine-tuned choice of parameters. Then the results are very 
similar to those for spinless fermions. For weak and intermediate impurity strengths 
the asymptotic low-energy behavior is approached only at rather low scales, which 
are accessible only for very large systems. This slow convergence observed already 
for spinless fermions holds also in the absence of two-particle backscattering. 

Interesting further extensions of the fRG for impurities in Luttinger liquids 
include the investigation of non-equilibrium phenomena and the analysis of disorder. 
For isolated impurities the influence of impurities on the vertex renormalization 
is irrelevant for the asymptotic low-energy or long-distance behavior, although it 
may contribute quantitatively at intermediate scales. For disordered systems with a 
finite impurity density the influen ce of the latter on the two-particle vertex is crucial 
and must be taken into account loiamarchil l2004| . In principle this is possible by 
computing the vertex flow with full propagators, which contain the renormalized 
impurity potential via the self-energy. A further challenging extension concerns 
the inclusion of inelastic processes. They appear at second order in the interaction 
and can be included in the flow equations by inserting the two-particle vertex into 
the flow equation for the self-energy without neglecting its frequency dependence. 
Finally, a flexible microscopic modeling feasible with the fRG approach allows for a 
more realistic description of contacts and leads in experimental systems. 
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A Evaluation of vertex flow for spin-| 
fermions 

A.l Functional RG 

Here we present a detailed derivation of the flow equations for the two-particle 
vertex T^. Starting from the flow equation ()4.5H) we insert on the right-hand side 
the parametrization ()4.47j) for and ()4.49j) for F^. The flow of the singlet vertex 
F^ , , is computed for the three choices of {k'l, ki, fcg) which yield the flow 
of 

fi's2' fi'sS' dti corresponding to Eqs. ()4.42|) and ()4.43p. while the flow of the triplet 
vertex F^ ^ , is evaluated only for (/c'^, /cg, /ci, /C2) = {kp, —kp, kp, —kp) as in 
()4.4Hl . which yields the flow of gf". For a = s2, s3, s4, t we obtain 

uj=±A ■ 

with 



(2P, + f/,-2^of/^cosp)2 , {2Ui + U,-fio{Ui + P,)cospy 



(4 - fiDiUj - PJ^ sin^p + 6(2 + /io cosp)U',i2Ui + U,) 

6(2 + /io cosp)fioU^{U', + PJ cosp - (4 - - PJ sin^p 

3(2cosp + /io)2?7'? , ((/i^-2)t/^ + f/, + 2(?7^ + Pjcosp)2 



+ 



6(/ig - 2 - 2 cosp)^;((/ig - 2)U'^ + U^ + 2{U'^ + PJ cosp) 
4(^^-e°-.,)(^^-epV,) 

3{^il-2-2cospfU'l 



(A.2) 
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A.l Functional RG 



(2^7; + - fipjU', + PJ COSP + 2(^7; - PJ sinpsinfc^)^ 

+ 2(zc. - 

6(2 + cosp - 2 sin p sin kF){2U'^ + U, - HoiU'^ + PJ cosp)f// 

12(2 + /io cosp — 2 sinpsin kp)U^{Ui. — P^) sinpsin /ci;' 
2(^^-e°)^ 

3(2 + yUocosp — 2sinpsin/ci?)^6'"'^ 

2(^^ - 

2(4 - /ig)^-^in^p (2[/; + - /io([/; + PJ cosp)^ 
(4 - ^o)iU', - Ps? sin^P - 2(2 + /ip cosp)(2y^ + UM 
2(/io(2 + /XQ cosp)(L/^ + PJ cosp - (4 - /ig)([/: - PJ sin^p)[/; 

5(/io + 2cosp)^ ((;,2_2)^/ + ^^ + 2([/^ + Pjcosp)^ 
4(^^ - O'U'l - - e°+.J 

2(/ig - 2 - 2 cosp)((/ig - 2)U'^ + U^ + 2{U'^ + PJ cosp)^/ 
4(^.;-L°. )(...- J 



(A.4) 



5(/ig-2-2cosp)2f/'? ^^^^^ 
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A Evaluation of vertex flow for spin-| fermions 



Here = —2 cos k — with /iq = —2 cos kp is the bare dispersion relation relative 
to the bare Fermi level. Since the functions fa{Pi ^) can be written as rational func- 
tions of cosp and sinp, the p -integral of Eq. (jA.lj) can be carried out analytically 
using the substitution z = e*^ and the residue theorem, in analogy to the spinless 
case described in Sec. 14.2.31 The resulting differential equations for the momentum 
space couplings read 



/2 



7(A) 



l\2 



7r(4 - (/io + «A)2) V iK 

(a^o + ikmU's + Us? + 2^^l{U', + P,f) + 4/xo(2^^ + u;){u'^ + PJ 

2(4-(/io + 2A)2) 

^g(/io + zA)=^([/^ + P,Y + 24(/xo + zA)(4[/^ + 2[/, - + P,))U', 

8(4-(/io + 2A)2) 

12/io([/, - 2P^)U[ - 3((/io + zA)(/ig + 4) - S^ji,)U'] 
2(4-(/io + ^A)2) 

(/io + 2A)(((/ig - 2)^- + + A{U', + PJ^) 

2/xo([/^ + Pj((/ig - 2)[/^ + U,) - 6/io((/ig - 2)P, - UM 
(/io + ihf - /io 

3(/io + iA)((/ig - 2)((/ig - 2)[/^ + [/J - 4([/: + PJ)[/; 
(/io + iA)2 - fil 

3((/io + 2A)((/ig-2)=^ + 4)-4/io (/ig - 2))t/' 

2((/io + 2A)2-/i2) 



/2 



(/iO + 2A)(f/'^ + P 



(/.o + ^A)(/.2(f/^ - PJ2 _ i2/.2p^f/; _ 48U'M + 24f/^ 



/2n 



(A.6) 



dgs3 
OA 



/2 



+ Re 



7(A) / 4f/ 



TT 



/2 



+ 



27r 



^J(^i + n + 6^D 
/io + ^A 
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A.l Functional RG 



4(/xo + iA) 



4 - (/io + lAy 



4(2P, - U^Y + /io(/io + tA){2U: - U,)i2Ui - U, - 12U^) 
4/io(4-(/io + zA)2) 



(A.7) 



dgs, _ 8U'i + (/ig - 2){U'i + - 3U't + QU'M) + AU'.P^ + 12P,y; 
OA 



+ Re 



7(A) 



47r 



2(/io + ^A)P,(t/: + 3t/;) 



+ 



7r(4 - (/io + zA)2) 



iifil - 2)P, + U, + ((/io + iA)/io - 4i7(A) sinA:^)f/^) 



/\2 



+ 4z7(A) sin kp - 4) 



/ioA2 

{^^o+^A)m:+u,r+^^m+p,,r) 



(Ato(Ato + ^A) 



+ 



+ 



4-(/io + ^A)2 

4/xo(2[/^ + ^J(^^ + PJ - 6(/io + tA)iAU'^ + 2[/, - /ig([/^ + PJ)^; 

4 - (/io + iA)2 

12/io(f/, - 2PJf/; - 3((//o + ^A)(4 + fil) - 8/io)f/1 



4-(/io + ^A)2 



(A.8) 



% (4-/ig)([/^P,-P,^;-2[/-g) 
9A Air 

, p r 7(A) / (/io + zA)(4-/ig)P,([/^-[/;) 

L 7r(4 - (/io + zA)2) I 2 

(4 - /ig)(4 - (/io + tAf)U", 2/io(2^^ + ^J([/^ + PJ 
^A (4 - (/io + ^A)2) 

(/io + ^A)((2[/- + [/j2 + /ig([/- + Pj2) 
2(4-(/io + 2A)2) 

((/io + tA)iAU'^ + 2[/, - filiU: + PJ) - 2/io(^, - 2PJ)^; 
4-(/io + ^A)2 
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A Evaluation of vertex flow for spin-| fermions 



5((/io + zA)(/ig + 4) - 8^o)f/1 (/io + ^m^^l " ^'s + U,f 



+ 



2(4-(/xo + zA)2) 2((/io + zA)2-^2) 

2(/io + zA)([/: + PJ^ + 2/io([/; + Pj((/ig - 2)^; + 
(/io + zA)2 - /ig 

(/XQ + ^A)((/xg - 2)((/ig - 2)U', + [/J - 4(^; + PJ)^; 
(/io + iAy - /io 

2/.o((/ig - 2)P, - UM 5(/io + ^A)((/i2 - 2)2 + 4)[/' 

(/io + ZA)2 - /i2 

10/io(/i2 - 2)[/^g 
(/io + iAy - fil 



/2 



2((/io + zA)2-/,2) 



(A.9) 



with 



7(A) 



'1- 



(/io + «A)^ 



(A.IO) 



Using the hnear equations ()4.48|1 and ()4.5()j) to replace ga by the renormahzed real 
space interactions on the left-hand side of the flow equations, we obtain a complete 
set of flow equations for the four renormalized interactions Uj-, U^, U!,, and P^ of the 
form (I13D). 



A. 2 One-loop g-ology calculation 

In the low-energy limit the flow of the two-particle vertex F^, as described in 
Sec. lA.ll reduces to the one-loop flow of the q- ology model, the general effective 
low-energy model for one-dimensional fermions js61vomlll979| . In the (/-ology ap- 
proach interaction processes are classified into backward scattering {gi±), forward 
scattering involving electrons from opposite Fermi points {g2±), from the same Fermi 
points ((74_l), and umklapp scattering ((/3_l). All further momentum dependences of 
the vertex are discarded. 

The g-o\ogy couplings are related to the momentum space couplings gs2, gss, gsA 
and gt by 

gi± = ^{gs2 - gt) 
g2± = l:{gs2 + gt) 
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A. 2 One-loop g-ology calculation 



9s3 

9.± = f (A.11) 
and to the real space couplings by 



9l± = ^{d^s2 - 


f 2r,3 + r,4 - (4 - 




92± = ^{f^l^s2 - 


f 2r,3 + r,4 + (4 - 




93± = ^(-4r,2 


■~ 2rs3 + Tgi) 




94± = ^(4rs2 - 


(2 - /i2)r,3 + r,4) 




respectively. 







(A.12) 



The flow equation for fj4.51|) reduces to the standard one-loop g-o\ogy cal- 
culation SolyomI [l979| . once the dependence on the internal momentum p on the 
right-hand side of the flow equation is neglected. Applying the above parametri- 
zation, we obtain a complete set of flow equations for gi± for i = 1, . . . , 4 of the 
form 

^ = 7^ [9i±92±PP{0) + 9i±94±PH{0) + 2g^^{g2± - gi±)PH{2kF) ] 

OA ZTT 

^ = ^ [li9l± + 9!±)PPiO) + 9A±i9i± - 92±)PHiO) + gl^PH{2kp) ] 
^ = ^ [93±9a^PP{^) + ^?3±(2^72± - 9iA)PH{n) ] 

^^'^ ^ [9l±PPi2kF) + hgl^ + 2g^^g2± - 2gl^ + gl^)PH{Q) ] , (A.13) 



dk 27r L^^^ ^ ' 2 

where 



1 



ip+q/2 — i-p+q/2 



PffM = ±Vr*^ (A.14) 

The initial conditions are g^l = g^l = U + (/ig - 2)f/' and g^l = g^l = U + 2U'. 
The integrals in Eq. ()A.14|) can be computed analytically using the residue theorem, 
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A Evaluation of vertex flow for spin-| fermions 



as discussed in Sec. lA.ll for g = and q = 2kp we obtain 

4 / Z7(A) 



PP(0) = --Re( 
PH{0) =4Re 



4 - (/io + iA) 



(/io + iA)7(A) 



PHi2k,) = 4 Re ( ^^ ^ ^) , (A.15) 



,(4-(/io + 2A)2)2, 

4 / 7(A) \ 

PP{2kF) - — Re _ ^ iA)2) (4 - /io(Aio + ^A) + 4i7(A) sin kp) ) 

(/io + 2A)7(A) 

,(4-(/io + 2A)2)(/i2-(^0 + ^A)2), 

where 7(A) is defined by Eq. (jA.lOj) . 

For /U 7^ the above equations reduce to 

PP(0) ^ 
PH{0) 
PP{2kF) 



Av/4^ 
4(2 + /ig)A 

8A 



P/J(2M = i (A.16) 

Av/4 - /i2 

in the hmit A — 0. The resulting flow equations in the low-energy limit are 
dA Avrv/I^^ 

%^ = — 4^ (A.ir) 

OA 2A7ry4^ ^ ' 

for fl'ix and g2±, whereas the flow for g3j_ and g4± vanishes. The solution reads 
, _ U+{f,l-2)U' 

gi± 



^2. - i + ^ + - ^0)^' j ' ^^-^^^ 

yielding the fixed-point couplings gi± = 0, g2± = -|- (6 — fil)U'), and g3± = 
U+{^Il-2)U',g^^ = U + 2U'. 
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B Bethe-ansatz calculation of for the 



p 



Hubbard model 



The Luttinger-liquid parameter Kp can be evaluat ed exactly from the Bethe ansatz 



solution for the one-dimensional Hubbard model |Lieb and WiJll968l |. The Bethe 
ansatz provides the energies of the ground and excited states as solution of specific 
integral equations. 

In the description of cr itical properties the "dressed charge matrix" is introduced 



Frahm and KorepinI ll99C| . This 2x2 matrix contains the effective renormalized 



coupling constants within and between the low-energy charge and spin sectors of the 
Hilbert space of the Hubbard model, and therefore directly determines all critical 
exponents. It is defined as 



(B.l^ 



In the absence of a magnetic field, Z is completely determined by its first element 




as 



i{k,) 

2 



(B.2) 



where fco is a cutoff determined by the particle density n Frahm and KorepinlllQQO | . 
^{k) obeys the integral equation 
4 r^o A 

^{k) = 1 + - / cos k'R{- (sin k - sin k') ) ^{k') dk' , (B.3) 



U 



fco 



with the kernel 



^ ^ 271 Jo l + ev ^ 
The cutoff momentum /cq is defined by 

fko 



p{k') dk' = n 



[B.4) 



(B.5) 



-fco 
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B Bethe-ansatz calculation of Kp for the Hubbard model 



where the integral equation for the ground-state charge distribution function p{k) 
reads 



14 r^o A 

p{k) = — + - cos A; / R{- (sin k - sin k') ) p{k') dk' (B.6) 
2vr U J_k^ U 



in the hmit L ^ oo . 

The compressibihty k is related to the dres sed charge matrix element ^(/cq) by 
'C^(^o) = 'n'VpK, where Vp is the charge velocity Frahm and KorepinI 199fl| . On the 
other hand, the Luttinger-liquid parameter Kp is given by Kp = 7TVpK/2. Hence Kp 
is determined by 

K,^^-^. (B7) 

Solving the integral equation for ^(/cq) with ko determined from equations ()B.5|) and 
()B.6|) . yields the exact Luttinger-liquid parameter. For the numerical computation 
of the kernel R{x) defined in Eq. ()B.4jl . the following expression is more convenient 



1 °° 9/ 

Similarly the chemical potential fi is determined. A comparison of the fRG 
results with the exact results allows a quantitative estimate of the effect due to the 
neglected frequency dependence of the two-particle vertex. The chemical potential 
can be derived from the ground-state energy via 

^=d^ = dkSdkJ ' ^^-'^ 

where the energy per lattice site is given by 

COS kp{k)dk. (B.IO) 

'ko 



Using equations ()B.5|) and ()B.6|) . the solution of the following integral equations 
de ^ , , ^ , dp{k) 



dk 



/Ko 
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with 



dpik) 4 4 4 

-rr, — = 77 p^h) cos A; [ i? ( — (sin - sin fco) ) + ^ ( 77 (sin k + sin ko) ) 1 
okQ U U U 



4 /■'^o 4 

+ — cos k ^{tt (^^^ ^ ~ ) 



(B.12) 



provides the exact chemical potentiaL 

Analogous expressions can be derived for the spinless fermion model 



Haldane 
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